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Abstract 

We analyse the structure of Yukawa couplings in local SU(5) F-theory models with 
enhancement. In this setting the symmetry is broken down to SU(5) by a 
7-brane conhguration described by T-branes, all the Yukawa couplings are generated 
in the vicinity of a point and only one family of quarks and leptons is massive at 
tree-level. The other two families obtain their masses when non-perturbative effects 
are taken into account, being hierarchically lighter than the third family. However, 
and contrary to previous results, we hnd that this hierarchy of fermion masses is not 
always appropriate to reproduce measured data. We hnd instead that different T- 
brane conhgurations breaking to SU(5) give rise to distinct hierarchical patterns 
for the holomorphic Yukawa couplings. Only some of these patterns allow to ht 
the observed fermion masses with reasonable local model parameter values, adding 
further constraints to the construction of F-theory GUTs. We consider an E^ model 
where such appropriate hierarchy is realised and compute its physical Yukawas, 
showing that realistic charged fermions masses can indeed be obtained in this case. 
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1 Introduction 


The proposal made in []J-|4| to build realistic 4d vacua by means of GUT constructions 
in F-theory has undoubtedly generated a wealth of activity in the past few years. While 
from all the classes of models considered in the string phenomenology literature [5]-[^ 
none is a priori preferred to achieve a fully realistic model of Particle Physics, F-theory 
vacua present a number of conceptual and technical advantages that has allowed to make 
substantial recent progress on this front. 

For instance, type IIB/F-theory flux compactihcations contain the class of vacua where 
at present moduli stabilisation is best understood [8 -10 , and where most of the statistical 


analysis of string vacua has been applied 11 . In addition, the construction of F-theory 


compactihcations relies on complex geometry, and so allows to implement powerful results 
in algebraic geometry. In practice, this translates into very useful techniques that can 
be used to construct explicit examples of F-theory models, as well as to gain a global 
perspective of the set of vacua as a whole. 

But perhaps the key ingredient that highlights F-theory GUTs as a promising avenue to 
construct realistic vacua is the localisation properties of 7-branes and their consequences. 
Indeed, the fact that 7-branes localise gauge and chiral degrees of freedom lets us formulate 
the construction of F-theory GUTs in a bottom-up fashion [^. This in turn permits to 
express the basic features that the gauge sector of realistic model should contain in terms 
of a small internal region where the GUT helds are localised, and to compute many 
quantities of physical interest in terms of such local data. 

A good example of the latter is the computation of Yukawa couplings in F-theory 


GUTs. Through a series of works 13 -26 it has been realised that to extract the flavour 


structure of an F-theory model one may implement an ultra-local approach and compute 
its Yukawas by analysing small regions of the four-cycle Fgut where the GUT degrees of 
freedom are localised. More precisely one hnds that, if gauge helds are localised in Fgut 
and chiral helds at complex matter curves Sj inside S'gut, then holomorphic Yukawas 
can be computed by looking at the points of intersection of such matter curves. Physical 
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Yukawas, on the other hand, can be computed ultra-locally if the internal wavefunctions 
for the chiral helds are sufficiently localised in a region near such point of intersection]^ 
An important outcome of the analysis of Yukawas in F-theory is that one may easily 
engineer GUT models where the Yukawa matrices are of rank one, by simply imposing 


a topological condition on the matter curves 22 . This will automatically give a mass 
hierarchy between one family and the rest. The masses of the two lightest families can 
then be generated when taking into account the effect of an Euclidean D3-brane instanton 
on a different four-cycle, along the lines of 20 Notice that rank one Yukawas are not 


exclusive of F-theory models (see e.g. 43,44 for type II examples). However, the F-theory 
framework does allow to compute them systematically for a wide class of models by means 
of the ultra-local approach. Quite remarkably, this remains true even after we include the 


non-perturbative corrections of 20 


In this spirit, the computation of Yukawa couplings in the presence of non-perturbative 


effects has been carried out in 24 26 . In particular, refs. 25 and 26 respectively 
analysed the fermion mass hierarchy developed for down-type 10 x 5 x 5 and up-type 
10 X 10 X 5 couplings in SU(5) models, which become MSSM Yukawas once that hyper¬ 
charge flux effects breaking SU{5) —)■ SU{3) x SU{2) x U(l)y are taken into account. 
In both cases it was found that a family hierarchy of the form (1, e, e^) is generated for 
fermion masses, with e a small parameter measuring the strength of the non-perturbative 
effect. Such hierarchy is already present at the level of the holomorphic Yukawas, and 
allows to £t empirical data upon taking e ~ 10“^ and including the dependence of physical 
Yukawas on the worldvolume fluxes threading Fgut- 


The computations carried out in 25 and are independent from each other because, 
in principe, down-type and up-type Yukawas can be generated at very different points of 


the GUT four-cycle S'gut- However, as pointed out in 45 experimentally the GKM matrix 
describes non-trivial correlations between U and D-quark Yukawas, and this strongly 
suggests that in a realistic model these two points should be in the same neighbourhood 
of S'gut, so that they experience the same local geometry and worldvolume flux densities. 

. For other applications 


^For a more precise statement in terms of the notion of local chirality see 
of seven-brane wavefunctions see 


27 


28 


33 


For different approaches to the generation of Yukawa hierarchies in F-theory GUTs see e.g. 


21034 


42 
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A model in which the two Yukawa points are very close to each other or even coincide 
is very attractive from the bottom-up perspective, as one is then able to compute the 
whole set of Yukawas with the mere knowledge of the local patch of S'gut- In this sense, 
even more appealing is the case where further matter curves intersect at a single point 
PEg G S'gut, such that the singularity of the F-theory elliptic hbre enhances to Eg at that 
point. One would then be able to compute ultra-locally the whole set of couplings which 
are most relevant for the GUT gauge theory by analysing a local patch around pEg , which 
is usually dubbed the point of Eg in F-theory (45) . 

Roughly speaking, the aim of this work is to apply the scenario of to F-theory 
models of SU(5) unihcation with a point of Eg where all Yukawa couplings are generated. 
More precisely, we consider F-theory models that can be locally described in terms of an 
Eg symmetry higgsed down to SU(5) by the 7-brane position background. This region 
of Eg symmetry contains both Yukawa points Pup and Pdown, which may be coincident or 
not, and the matter curves intersecting at them are chosen in such a way that we obtain 
rank one Yukawas at tree level. The question we address is then if, by taking into account 
non-perturbative effects, a realistic hierarchical pattern of Yukawa couplings is generated 


such that it allows to £t experimental data, in the same spirit of 24-26 


As mentioned, the hierarchy that has so far allowed to £t empirical data is of the 
form (l,e,e^), at least for reasonable values of e and local model parameters like flux 
densities. Therefore we take this hierarchical pattern as the guiding principle to achieve 
a realistic fermion mass spectrum in the context at hand. Remarkably, this hierarchical 
structure can already be seen at the level of holomorphic Yukawas, which depend on very 
few parameters of the local model. As a result, the above criterion is very robust and can 
be applied even without specifying the oftentimes complicated 7-brane flux background. 

As we will see, while one may construct several Eg models with rank one Yukawas at 
tree level, the addition of non-perturbative effects does not always yield a hierarchy of the 
form (1, e, e^). In fact, from all the models that we have analysed only one choice generates 
the desired hierarchy for both types of Yukawas and at the same time contains interesting 
mechanisms for realistic p-term and neutrino masses. Of course having this hierarchy 
at the holomorphic level does not guarantee that one can reproduce the whole set of 
empirical data related to Yukawa couplings. Hence, once selected the most promising 
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Eg model we proceed to describe it in full detail and to compute its physical Yukawa 


couplings. We find that, similarly to the results in 24 26 , within the MSSM scheme 
fermion masses at the GUT scale can be £t for e ~ 10“^ and large values of tan f3. This 
result is valid for the families of quark and leptons over which we have good control given 
the approximations taken in our analysis, namely the two heaviest families. This not only 
applies to the fermion mass spectrum but also to the quark mixing angles. We determine 
the latter in terms of the separation of the two Yukawa points Pup and Pdown, providing 
precise formulas that illustrate previous statements in the literature. 

The paper is organised as follows. In section we review the construction of local 
F-theory GUTs and how hierarchies of Yukawa couplings arise via non-perturbative ef¬ 
fects. In section we analyse several local Eg models and compute their Yukawas at 
the holomorphic level, selecting one model based on the criterion described above. In 
section]^ we describe this particular model in detail and in section]^ we solve for its chiral 
zero modes wavefunctions and compute the normalisation factors that render their kinetic 
terms canonical. Such normalisations are the missing ingredient to compute the physical 
Yukawas, whose hierarchies are analysed in section We conclude in section 

Several technical details have been relegated to the appendices. Appendix [X| contains 
the details and notation regarding the Eg Lie algebra used throughout the rest of the 
paper. Appendix contains the details of the computation of holomorphic Yukawas for 
the models of section Appendix discusses the notion of local chirality applied to the 
Eg model of section while appendix spells out the computation of its holomorphic 
Yukawas, and appendix its zero mode wavefunctions in a real gauge. 


2 Yukawa hierarchies in F-theory GUTs 

The standard description of F-theory GUT models [T]-[^ (see for reviews) is done in 
terms of a Galabi-Yau fourfold which is elliptically hbered over a three-fold base B. The 
hbration is such that the fibre degenerates over a 4-cycle S'gut, with a hbre singularity 
whose Dynkin diagram corresponds to Gqut- At certain 2-cycles S within 5gut the 
fibre may display a higher singularity type, signalling the presence of chiral multiplets 
charged under Ggut and localised at such matter curve. The precise 4d chiral matter 
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content of the model depends on the fonr-form flnx G 4 threading Sgut which, if chosen 
appropriately, can break Ggut to the snbgronp 517(3) x SU{2) x 17(l)y. Finally, the 
conplings among different chiral mnltiplets will depend on the intersection pattern of the 
corresponding matter cnrves. Hence jnst by knowing the local geometry aronnd 5 gut one 
may see which Yukawa couplings may be generated in the effective 4d GUT theory. 

While such a picture is rather compelling it is not the appropriate one for the actual 
computation of Yukawa couplings. Instead, it proves more useful to work with an alter¬ 
native description of the degrees of freedom localised at 5gut- Namely, one can use a 
8d action related to the 7-branes wrapping 5 gut and those intersecting them. Such an 
action is defined on a 4-cycle S and in terms of a non-Abelian symmetry group G D Ggut, 
and is such that the 4d gauge theory described above can be obtained upon dimensional 
reduction. In particular, Yukawa couplings can be obtained from the following superpo¬ 
tential 


W = mt / Tr (F A 4)) 


( 2 . 1 ) 


where m* is the F-theory characteristic scale, F = dA — iA A A is the field strength of the 
7-branes gauge boson A, and $ is the so-called Higgs field: a (2,0)-form on the 4-cycle S 
describing the 7-branes transverse geometrical deformations. Both A and 4) transform in 
the adjoint of the initial gauge group G, which is nevertheless broken to a subgroup due 
to their non-trivial profile. In particular, the prohle (4)) is such that it only commutes 
with the generators of Ggut in the bulk of ^gut, while on top of the matter curves of 
5 gut it also commutes with further roots of G. The background profiles (4)) and (A) 


cannot be arbitrary, but they must solve for the equations of motion that arise from (2.1) 
and the D-term 


D = [ a;AF + ^[4),4>t] 

JQ ^ 


( 2 . 2 ) 


where uj stands for the fundamental form of S. Given the background profiles for 4> and 


A one can use (2.1) and (2.2) to solve for the internal wavefunction of their zero mode 
fluctuations representing the 4d matter fields. Finally, these wavefunctions can be plugged 


back into (2.1) to compute the precise value of the Yukawa couplings of the model. 

In this setting, two important results that apply to the holomorphic part of the Yukawa 
couplings are i) they do not depend on the prohle of the worldvolume hux F and ii) they 
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only depend on the local geometry around the point p where the involved matter curves 


intersect 17 . Thanks to this, in order to compute holomorphic Yukawa couplings it 
suffices to describe the prohle (d*) around a neighbourhood Up C Sgut of the intersection 
point p, and take G = Gp to be a symmetry group just large enough to contain Gqut 
and describe the matter curves intersecting at p. This ultra-local approach is also valid 
to compute physical Yukawa couplings if one assumes that the wavefunction profile for 
the corresponding 4d chiral fields is peaked within Up, something usually achievable due 
to the localisation properties of matter curves and of the worldvolume flux F threading 
them. 

This approach to compute Yukawa couplings has been mostly developed in the case 
where Gqut = SU{5). There we have that the up-type Yukawa couplings 10 x 10 x 5 
can be described by taking Gp = Eq or larger, while the down-type couplings 10 x 5 x 5 
require at least of Gp = SO{12). While in principle one may consider several intersection 


points of each kind, it was proposed in 13 a scenario where all up-type Yukawas are 
generated from a single Yukawa point pup, and all down-type Yukawas from pdown, on 
the grounds that it is then natural to obtain that one fermion family much heavier than 
the other two. In fact, what one finds for the down-type couplings is that the Yukawa 
matrix has rank one, and hence only one family of quarks and leptons becomes massive 
[It] . This one-rank result is rather robust in the sense that it only depends on the 
topological intersection pattern of matter curves, and so deforming the divisor S'qut or its 
worldvolume flux F will not change it. Nevertheless, non-perturbative effects associated 
to other 4-cycles Snp C B will affect the result and increase the Yukawa rank from one to 


three 20 


Indeed following the discussion in 20 , if Snp hosts a 7-brane with a gaugino condensate 
or a 3-brane instanton with the appropriate number of zero modes then the S'qut 7-brane 


superpotential (2.1) will be modified, obtaining 


W = mt / Tr (F A 4)) + e ^Tr (F A F) (2.3) 

Js 2 

where e measures the strength of the non-perturbative effect, and do is a function that 
depends on the embedding of the 4-cycle S'np (we have that 6*o = (47r^m*)“^[log h/ho]z=o, 
with h the divisor function such that Snp = {h = 0} and ho = Jgh). Further corrections 
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that depend on the derivatives of h normal to S do also appear, as well as corrections 
at higher powers of e, see for explicit expressions. Nevertheless, for realistic 

models the increase of the Yukawa rank and the generation of hierarchies can already be 


seen from the leading correction shown in (2.3), with which we will work in the following, 
while the remaining contributions are rather suppressed. 

The analysis of down-type Yukawa couplings for Ggut = SU{5), Gp = SO{12) and 


with the superpotential (2.3) was carried out in ref. 25 . It was obtained a fermion 
mass hierarchy of the form {0{l),0{e),0{e‘^)) that is already present at the level of 
the holomorphic Yukawas. Such robust hierarchy allows to £t the experimental values of 
quark and lepton mass ratios once run to the unification scale. For this one needs e ~ 10“'^ 
and to take into account the worldvolume flux dependence of the physical Yukawas. In 
particular, their dependence on the hypercharge flux FV allows to obtain a realistic ratio 
for the r and 6-quark Yukawas and for the different mass quotients between the second 
and third families of D-quarks and leptons. Being allowed to fit all these data is a clear 
improvement with respect to classical 4d held theory models of SU(5) unification. 

In principle one can implement the same strategy to obtain a realistic spectrum of 
up-type Yukawa couplings. However, because the coupling 10 x 10 x 5 involves two 
fields with the same quantum numbers one cannot achieve a tree-level rank-one result for 
these Yukawas with the simple intersection of three matter curves. Instead, one needs to 
consider more involved matter curve geometries which, in terms of the 7-brane position 
field <F can be described by a non-Abelian background profile for (<h), usually dubbed 


T-brane 19,^ (see also 46-48]). 

The presence of T-branes does complicate the 7-brane wavefunction equations but, as 


shown in 26 , one can still analyse the case of up-type Yukawas by taking Ggut = *S'f/(5), 


Gp = Eq and the superpotential (2.3). It was found in this reference that the inclusion of 
non-perturbative effects also modifies the tree-level rank result to a hierarchical structure 
of the form (0(1), 0{e), 0(e^)) for the holomorphic Yukawa eigenvalues. Finally, realistic 
values for the top-quark and for the quotients of t/-quarks can also be achieved by again 
taking e ~ 10“^ and values for the worldvolume flux densities very similar to those needed 
around Pdown to £t the down-type Yukawa data. 


These previous results suggest that one may naturally obtain a realistic mass spectrum 











by considering pup and pdown in the same neighbonrhood of Sgut, as this wonld explain 
why local parameters like flnx densities are similar aronnd both points. In practice, this 
amonnts to consider a symmetry gronp large enongh to describe the whole set of matter 
cnrves containing the MSSM chiral content, which selects either Gp = Ej or Gp = Eg. 


In fact, it was proposed in 45 that considering pup and pdown to coincide in a point of 
Eg enhancement would account for all the flavour hierarchies observable in the Standard 
Model. In this paper we would like to investigate if that is indeed possible in the scheme 
discussed above, namely where the hierarchies are obtained from perturbing a tree-level 
rank-one result via non-perturbative effects, as encoded by the corrected superpotential 

@0 


The set of Eg models that can accommodate the matter spectrum of the MSSM is 
rather rich, and it is usually classified in terms of the matter curves present in S'gut- 
This classification is particularly powerful whenever the matter curves content admits 
a spectral cover description, which in particular means that there is an underlying Eg 


structure globally dehned over S'gut 49 - 53 or even through the whole threefold base 54 


In our ultra-local approach such spectral cover description is not necessary, in the sense 
that the Yukawa couplings will depend on the profiles of $ and E near the Yukawa point 
P = Pup = Pdown- In particular, the holomorphic Yukawas will only depend on the profile 
of <I> around p. As the hierarchy of couplings is already captured at the holomorphic level, 
we will proceed to classify our Eg models based on their local profile (d*). Finally, just 
as in the Eq case we need the presence of T-branes in order to have a tree-level rank-one 


up-type Yukawa matrix, so the models to consider will be T-brane prohles for the held <I> 
around the Yukawa point p of Eg enhancement]^ 

^To be precise, our scheme only assumes that the hierarchies present in the up-type and down-type 
Yukawa matrices are all due to a rank-one tree-level superpotential corrected by non-perturbative effects, 
while the hierarchies in the neutrino sector could be due to a different mechanism. In this sense it is 
equally interesting to explore the presence of hierarchies in points of ify enhancement, which will be 
discussed in a separate publication. 

^Strictly speaking due to the T-brane profile (<i>) does not vanish and so the symmetry is not enhanced 
to Eg at any point of the local model. In particular at the point p where all matter curves meet and all 
Yukawas are generated we will get a fibre singularity enhancement but we do not expect to recover a full 
Eg singularity. For the sake of simplicity, we will abuse of language and still refer to this p as a point of 
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To summarise, the findings of 25,26 suggest one can describe the flavour hierarchies 


of the Standard Model via an F-theory SU(5) local model with a point where all the 
Yukawas originate from. The geometry near such point should be described by an T- 
brane profile for $, and non-perturbative effects should be taken into account to increase 
the rank of the Yukawa matrix from one to three. Since the holomorphic Yukawas already 
detect the fermion mass hierarchies, one can already see at this level whether a specific 
T-brane model is promising for reproducing empirical data. In the following section we 
will consider different classes of T-brane backgrounds and analyse if they give rise to the 
appropriate flavour hierarchy of fermion masses. 


3 SU{5) models with Eq enhancement 

In this section we will present a set of local SU(5) F-theory models that can be described as 
an Eg theory higgsed by a T-brane background. Each of these models has the appropriate 
structure of matter curves so that they can embed the full content of the MSSM chiral 
spectrum, with only one massive family at tree-level. The remaining families of quark and 
leptons will become massive due to non-perturbative corrections, but then we find that 
one may get a hierarchy of masses either of the form (l,e, e^) or of the form (l,e^,e^). 
Since e is a very small number that measures the strength of a non-perturbative effect, the 
latter hierarchical pattern is very unlikely to reproduce empirical data, while the former 


has already been shown to be adequate in simpler local F-theory models 25,26 


As discussed in the previous section, to discover the hierarchical pattern that non- 
perturbative effects give rise to it suffices to compute the holomorphic piece of the Yukawa 
couplings. This will greatly simplify the analysis of this section, as these couplings can be 


computed via a residue formula 17,22,25,26 . Finally, for the sake of clarity we will only 


display the basic features of each model, their structure of matter curves and the final 
result for the holomorphic Yukawas, leaving the computational details to appendix 


Eg, in the model. 
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3.1 T-branes and matter curves 


One crucial feature of an F-theory local model with respect to the computation of Yukawa 
couplings is the profile for the 7-brane Higgs held <F in the vicinity of the Yukawa point. 
As mentioned above, obtaining a third family much heavier than the other two naturally 
selects a T-brane prohle for <F, which then specihes an appropriate local structure of matter 
curves. In the following we will briehy discuss the relation between T-brane prohles and 
matter curves, as they will be used when discussing each model. For a more thorough 


discussion on this subject we refer the reader to 22,26 


A T-brane background is specihed by a particular conhguration for the Higgs back¬ 
ground (<F) that does not commute with its adjoint, namely [(*h), (‘h)l] 7 ^ 0. In this class 
of backgrounds the identihcation of the matter curves can be subtle because, unlike in the 
case of commuting Higgs field, there will not be an enhancement of the symmetry group 
in a complex curve within S. Nevertheless, it is still true that some additional roots of 
the algebra of Gp will commute at specihc complex codimension one loci, and this allows 
us to identify the matter curves as these particular loci. 

In order to detect the structure of matter curves it proves useful to work with matrix 
representations of the Higgs field in the algebra dehned such that 0gut © 0 _l is a 
maximal subalgebra of Qp = Lie(Gp). Let us consider the case of interest in this paper, 
namely 0 p = eg and 0gut = -SUs . Then we have the well-known maximal decomposition 


,GUT 


Cg D SU5"" ©SUg 

248 ^ (24,1)® (1,24) ©((10,5)© c.c.)© ((5,10)© c.c.) 


(3.1) 


Since the Higgs profile (<h) belongs to the adjoint of Cg and by construction commutes 
with it will only act non-trivially on the each of the representations TZ of q± = su^ 


that appear in (3.1). This action can be expressed in terms of a matrix $ 7 ^ such that 


[(<F),77] = ^tiTZ (see e.g. 26 , section 3) so whenever the determinant of $ 7 ^ vanishes an 
element of TZ will be commuting with (<h). Finally, given that {TZq\jt,TZ) C 248, this will 
indicate that we have a zero mode transforming as TT-gut in the locus where det $ 7 ^ = 0 , 
and so the corresponding matter curve. 

Interestingly, these facts allow to express the structure of matter curves in terms of 
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the spectral surface of the Higgs held, which is dehned a^ 


y, z) = det(<h 7 e - zI)=Q 


(3.2) 


for each of the matrices $ 7 ^ associated to (<h). Following [^, we say that $ 7 ^ is recon- 
structible if its spectral surface is a non singular algebraic variety, and that it is block 
reconstructible if it has the structure of a block diagonal matrix such that every block is 
reconstructible. As the property of reconstructibility is independent of the representation 
TZ we then say that the Higgs held is block reconstructible, and in this case the whole 
information of {$) is carried by its spectral surfaces]^ 

Now it is easy to see how the pattern of matter curves can be encoded in the spectral 


surface (3.2): when the Higgs held is block reconstructible its spectral surfaces will be the 


product of polynomials whose zero locus is a non-singular algebraic variety, and there will 
be a one to one correspondence between these varieties and the matter curves in a specihc 
representation. Hence, the presence of several matter curves will induce a splitting of 
the spectral surface into irreducible polynomials, the number of factors of this splitting 
matching with the number of matter curves. 

In the following we will present a number of local Fg models whose local spectrum 
of matter curves can be detected by means of the above considerations. For the sake 
of simplicity, we will focus on models in which the Higgs held background (<h) is block 
reconstructible, since then we can classify our models by the number of matter curves 
near the Yukawa point. It would however be interesting to extend our set of examples to 
more general, non-reconstructible backgrounds. 


3.2 Catalogue of models 

We now proceed to describe several kinds of local Fg models with only one massive family 
at tree level. Such models are candidates to yield a realistic hierarchical fermion mass 
pattern after non-perturbative ehects have been taken into account although, as already 
advertised, this will not always be the case. To hnd out we will compute the holomorphic 


Yukawa couplings, which depend on the prohle for <F in the holomorphic gauge 17 . For 


’The following expression for the spectral surface holds if $ takes values in a u„ subalgebra of Qp. 


®As shown in 22 , SU{k) reconstructible T-branes correspond to spectral covers with monodromy 
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this purpose we only need to specify ($) as a linear combination of holomorphic functions 
multiplying the roots, following the notation of appendix As explained above we 
may also describe this background as a matrix $5 acting on the representation 5 of su^, 


which allows to find the local set of 10 matter curves via eq.(3.2). Since we are considering 


reconstructible backgrounds, the 5x5 matrix $5 will be automatically block diagonal, so 
we can classify our local models by the different dimension of each of these blocks. 

For simplicity we will only provide the basic data for each of the models that we 
have studied, leaving the details for appendix That is, we will describe the local set of 
matter curves that arise from the prohle for <F and the different possible assignments of 
the MSSM helds within them. Recall that in the absence of hypercharge flux the matter 
spectrum is organised in SU(5) multiplets, so for the purpose of computing holomorphic 
Yukawas we can consider that SU(5) is unbroken. Then to achieve rank one Yukawas at 
tree-level we need to have three copies of the matter representation 10 m within the same 
10-curve and three copies of 5 m in the same 5-curve. Finally the Higgs multiplets 5^/ 
and bo should be in 5-curves different from the one of 5 m and such that the couplings 
10m X 10m X 5(7 and 10m x 5m x br, are allowed. 

For each assignment we will present the structure of holomorphic Yukawa couplings 


that arise from the superpotential (2.3) with 


Oq = i{x9x + y9y) 


(3.3) 


and (x, y) parametrising the complex coordinates of the 4-cycle S. We will then discuss 
whether such structure accommodates favourable hierarchies to fit empirical data. In the 
next section we will provide a more detailed description of one of the models with such 
favourable structure, providing all the details that allow to compute its physical Yukawas. 


4-|-l models 

Let us first consider a holomorphic background for <F = ^xydx A dy of the form 

(‘f’xj/) = 2H2 -|- 3.H3 -|- 41^4) -|- 7 ji[E^ E2 mxE^ ) ( 3 - 4 ) 

where the notation and definitions that are used for the E^ roots are given in appendix [A} 
Here A = y^{hx—y) is a holomorphic linear function of (x, y) vanishing at the origin, which 
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is where the Yukawa point p will be located. By acting on the fundamental representation 
of su^ we obtain the matrix representation 




A m 0 0 0 

0 A m 0 0 

0 0 A m 0 




V 


m‘^x 0 0 A 

0 0 0 0 


0 

-4A 


(3.5) 


which displays a 4 + 1 block structure. The various matter representations and their 
matter curves are then the following ones 


- 10 sector 


5 sector 


10a : A"^ = m^x , 10b : A = 0 , 


5a : (3A)"^ = m^x , 5b : X^{m^x + 4A"^) = 0 , 


and it is easy to see that all the curves meet at the origin. 


This local model has already been considered in 23 , where it was found a rank one 


structure for the holomorphic Yukawas by using the tree-level superpotential (2.1). In 


the following we would like to extend this result by considering the superpotential (2.3) 


corrected by non-perturbative effects and providing the resulting holomorphic Yukawas 
up to order 0 {e^). 

As pointed out in (see also appendix in order to generate an up-type Yukawa 
coupling 10m X 10m x 5u it is necessary to assign the representation 10m to the 10 a 
curve and to the curve 5b. Because there are only two 5-curves, we will also consider 
that 5 d is also localised in 5b while the three copies of the 5 m are in 5a. 

With this setup one can compute the Yukawa matrices via a residue calculation. 
Schematically we hnd that 


/ 


Yu = 


0 

0 




0 eyi3 

ey22 0 

y £2/31 0 1/33 y 




(3.6) 
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/ 


^D/L — 


0 0 e 1/13 

0 e 1/22 0 


+ 0{^) , 


(3.7) 


y ^2/31 0 1/33 y 

where are order one numbers (detailed expressions for yij are found in appendix [B|). 
We then reproduce the results of in the limit e —)■ 0, while we see that for e 7 ^ 0 the 
rank of both matrices is increased to three. Finally, both matrices will have a hierarchy 
of eigenvalues of the form (0(1), 0(e), 0(e^)) so this model has a Yukawa structure which 
is favourable to reproduce the empirical data. 

Despite this favourable hierarchy, this model has the less attractive feature of having 
both up and down Higgses 5[/, in the same curve. Hence some particular mechanism 
should be invoked to prevent a large /i-term to be generated. Because of this potential 
drawback we will not consider this model in the following. 


3+2 models 

We next consider a Higgs background of the form 


i^xy) — —X + 2 .H 3 + 

^ +m{Ef + E'^+myE^ ) + m{EfQ+mxE^Q). 

(3.8) 

where again A = y‘^{bx — y). Its action 

on the fundamental of su^ is given by 



1 

coito 

m 

0 

0 

0 ^ 




0 


m 

0 

0 



$5 = 

fh?y 

0 


0 

0 

1 

(3.9) 


0 

0 

0 

A 

m 





0 

0 

mX'x 

>' j 



showing a 3 + 2 block structure. The various 

matter representations and curves are 

now 

- 10 sector 









g 

10 a ; —TpzA^ + m^y = 0 , lOj, : — vr?x = 0 , 


- 5 sector 


64. 


5a : m^y + = 0 , 5^ : rr^x^ = 

^ i 

5, ; A = 0. 




+ nrx 


2Xm^y - + 27m^2/)^ , 
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In this class of models we can assign the three copies of 10 m to either the 10a or the lOf, 
sector. If we assign the 10 m to the 10a sector then we need to assign to the 5a sector 
in order to have 10 m x 10 m x 5{/ Ynkawas which are singlets nnder S'f/(5)j_. Nevertheless, 
an explicit compntation shows that the holomorphic np-type Ynkawa conplings vanish for 


this arrangement. This vanishing resnlt is analogous to the one found in 23 for the E'j 
model studied in there, with the matter curves involved in the up-type Yukawas having 
a similar structure. We therefore see that this assignment of chiral matter to curves does 
not yield realistic Yukawas. 

The other possibility in this model is to assign 10m to the lOj, sector, which requires 
that 5f/ corresponds to the 5c sector. In addition one has to choose how to assign the 
representations 5 m and to the sectors 5a and 5b, having two possibilities. The final 
structure of the Yukawa matrices does however not depend on this choice. In both cases 
we find that 

^ 0 0 £2/13 ^ 

Yu = 0 £2/22 6 2/23 +C>(£^), (3.10) 

61/31 6 2/32 1/33 j 

^ 0 0 0 ^ 

Ydil= 0 0 £2/23 + C>(£^), (3-11) 

0 £2/32 1/33 y 

where again 2/ij are order one numbers whose explicit expression is given in appendix 
We see therefore that for this class of models the down-type Yukawa matrix does not have 
a favourable hierarchical structure. 


2-|-2-|-l models 

We finally consider a Higgs background of the form 
{^xy) = \i{Hi+2H2-2Hi)-\2{H2, + 2H4)+m{E^ +mxE^)+m{E2 +myE2) (3.12) 
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whose action on the fundamental of sug is 

/ 


m^y 


V 


0 

0 

fh 

—2Ai — A 2 
0 


0 

0 

0 

0 




(3.13) 


2(Ai + A2) J 


\i m 0 

m'^x Ai 0 

0 0 -2Ai - A 2 

0 0 

0 0 0 

and where now Ai and A 2 are two different polynomials of x, y which we shall take as 
Ai = fi^bx — y) and A 2 = yli^x — y). As we will see in the next section, taking Ai 7 ^ A 2 
with a slightly more general Ansatz will allow us to separate the two Yukawa points pup 
and Pdown from each other, introducing an interesting source of family mixing. 

The matter representations and matter curves are in this case 

- 10 sector 


10a : A^ — m^x = 0 , lOfe : (2Ai + A2)^ — m^y = 0 , 10c : Ai + A 2 = 0 , 


- 5 sector 

5a;Ai=0, 5b;2Ai + A 2 = 0, 5c : (3Ai + 2A2)^ - = 0 , 

5d : A 2 — fh^y = 0 , 5e : (Ai + A 2 )"^ — 2(Ai + \2)‘^{m^x + m^y) + {m^x — m^yY = 0 . 


so the amount of matter curves increases considerably with respect to previous models. 

In this case we can assign the representation 10 m to either the 10a or the 10;, sectors. 
Since both choices end up leading to the same results we will choose the first option, 
which hxes the 5c/ representation within the 5a sector. The up-type Yukawas then have 
the following structure 


/ 


Yu = 


0 

0 


0 epi3 

ei/22 0 




0{e^ 


(3.14) 


\ epsi 0 P33 j 

We then hud an eigenvalue hierarchy of the form (0(1), 0(e),0(e^)) and therefore a 
suitable hierarchical structure to £t empirical data. 

There are some possibilities now on how to associate the representations 5 m and 5// 
to the remaining matter curves, and this choice affects the down-type Yukawa matrix. 
We list here the possible choices and the resulting Yukawa matrices: 
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- Either 5 m is associated to 5^ and 5^5 is associated to 5e or the other way round. In 
the hrst case we hnd that the down-type Yukawa matrix has the form 

^ 0 0 0 ^ 

ko/L = 0 0 e?/23 +0(e^) (3.15) 

0 e?/32 2/33 y 

whose eigenvalues are (0(1), O(e^), 0(e^)). Therefore this assignment for the mat¬ 
ter helds does not lead to a good hierarchical structure for the down-type Yukawas. 
On the other hand, if we identify the 5 m with 5e, then it is not clear how to perform 
the analysis due to the fact that the matter curve is singular at the Yukawa point. 


- Either 5 m is associated to 5^ and 5^) to 5c or the other way round. In both cases 
we hnd that the down-type Yukawa matrix has the structure 


/ 


Yd/l — 


\ 


+ 0{e^). 


0 0 e 2/13 

0 e 2/22 0 

y ^2/31 0 2/33 y 

that has the favourable eigenvalue hierarchy (0(1), 0(e),0(e^)). 


(3.16) 


We then see that in the present 2 + 2 + 1 model there are two particular assignments 
of matter helds that yield a promising hierarchical structure for both up and down-type 
Yukawas. However, as we discuss next, one of the two has a more attractive structure for 
the p-term and neutrino masses, namely the choice of assigning 5 m to 5b and 5d to 5c. 


3.3 Comments on /i-term and neutrino masses 

One of the most attractive features of the Eg models is that it is possible to describe the 
masses of the neutrinos and the /i-term for the MSSM Higgs sector at the same time as 
the Yukawa couplings [^. Although in the following sections will not deal with them, 
let us briehy analyse here the structure of neutrino masses and /i-term in the case of the 
two 2 + 2 + 1 models with good hierarchical structures for the Yukawa matrices, in order 
to select one of them. 

In the case in which we assign the 5 m to 5c and the 5d to 5^ we hnd that there is 
a singlet under SU{5)gut that can give a coupling of the form 1 x 5;/ x 5 m and, after 
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breaking SU{5)gut down to the standard model gange gronp this will imply the presence 
of the following conpling in the snperpotential 


WdXH^LS, (3.17) 

where we called the singlet S and L the lepton donblet snperheld. This conpling, as 


analysed in 57 , corresponds to a Dirac mass for the nentrinos if we identify the singlet S 
with the right handed nentrino Nj^. With this assignment of matter cnrves however it is 
not possible to have a renormalisable /i-term for the Higgs fields. It is possible to generate 
a non-renormalisable /i-term nonetheless if we consider the interactions of the Higgs fields 
with modes coming from other matter curves. In particular when the fields in the 5e come 
in vector-like pairs the following couplings will be allowed in the snperpotential 

Wd XiHj<j) + \2HdS(j)^ + Acl><j)\ (3.18) 

where we called 0 any field in the 5e sector and (f)^ its conjugate, A is a mass term for 0 
and S' is a singlet. After integrating out 0 and (f)^ using their F-term equations we hnd in 
the snperpotential the following term 

A1A2 52 


W D 


A 




(3.19) 


which becomes an effective /r-term for the Higgs helds if the singlet S gets a non-vanishing 
vev. Note that this kind of non-renormalisable effective /i-term has already been consid¬ 


ered in 58 and can provide a solution to the /i-problem in the MSSM. 

In the second case, namely in the case we assign the 5m to 5^ and the to 5c, we 
hnd that it is possible to have the following coupling in the snperpotential 

W D SH^Hd (3.20) 

which becomes an effective /i-term if the singlet S gets a non-vanishing vev. This class of 
effective /i-term is particularly interesting because it can provide a mechanism for solving 


the /i-term problem in the MSSM 59 . However for this assignment of matter curves we 
hnd the feature that no masses for the neutrinos are possible if they are localised at the 
intersection of two 7-branes. Since one of the major motivations for studying Yukawa 
couplings at the point of is the generations of all couplings in the MSSM, including 
/i-term and masses for neutrinos, we will henceforth focus our attention on the hrst model 
discussed in this subsection and start analysing it in detail in the next section. 
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Summary 


To sum up, from all the models discussed in this section, we have found that the 2 + 2 + 1 


model specihed by (3.12) is the most interesting phenomenologically, in the sense that it 
yields a hierarchical structure of Yukawa couplings of the form (C>(1), 0{e), C>(e^)) for a 
specihc assignment of SU{5) representations to matter curves. Such assignment is 


10m 10a 5 m 5c 5(7 —>■ 5a 5() (3-21) 

which also exhibits interesting mechanisms to generate realistic neutrino masses and +- 
term. In the next sections we will analyse this model in detail, specifying a background 
that includes the appropriate worldvolume fluxes on each of the above matter curves 
(section]^. We will then compute the wavefunctions (section]^ and the physical Yukawas 
(section for this local Eg model, showing how empirical fermion masses and mixings 
can be £t upon an appropriate choice of parameters. 


4 An Eg model with hierarchical Yukawas 

Let us now consider in some detail the most promising of the Eg local models discussed 
above. In particular we will describe the 7-brane background for the last 2+2+1 model 
and show that it indeed incorporates a realistic hierarchy of Yukawa couplings. 

Describing the background of a 7-brane local model with Eg symmetry group entails 
specifying a Higgs held <h and gauge connection A, both valued in the algebra su^. If we 
want to preserve supersymmetry at the GUT scale, they must obey the F-term equations 


Oa^ = 0 (4.1a) 

F(°’2)=0 (4.1b) 


that arise from minimising the superpotential (|2.1|). Also, the D-term (2.2) should vanish 

(4.2) 


w A F + ^ 0. 


In order to hnd a solution to the above, one usually exploits the fact that the F-terms 
are invariant under the complexified gauge group, as opposed to the D-term which is 
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only invariant under the real group. More explicitly, any holomorphic Higgs together 
with = 0 automatically satishes the F-terms. This is referred to as a solution in 


holomorphic gauge 16,17 which is not physical since it does not obey the full set of 
equations of motion and the gauge held is not real. However, one can still extract useful 
information at the holomorphic level, such as the structure of matter curves and rank 
of the Yukawa couplings by solving a relatively simple algebraic problem. Finally, by 
performing a complex gauge transformation, we can bring the helds in a real gauge that 
satisfy the D-term. This last step is the most challenging as it requires solving a set of 
partial differential equations that become particularly complicated in models including 
T-branes. However, it is unavoidable if we want to obtain the kinetic terms and hence 
the magnitude of the Yukawa couplings. 

Following the above approach, we will hrst introduce the background for the Higgs held 
in holomorphic gauge and discuss the structure of matter curves and their intersections. 
We will then consider the background in a real gauge by imposing the D-term equation, 
which forces the introduction of non-primitive huxes. We will complete the description of 
the background by introducing additional huxes to achieve chirality for the MSSM helds as 
well as GUT symmetry breaking. This choice of background will yield an MSSM spectrum 


whose holomorphic Yukawa couplings are of the form (3.14) and (3.16), as we show by 
giving explicit expressions computed by means of a residue formula. The computation of 
physical Yukawa couplings and mixing angles is left for sections and 


4.1 Higgs background 

Holomorphic gauge 

The hrst ingredient necessary to dehne our local model model is the background Higgs 
held (<h) = {^xy) dx A dy that triggers the breaking of to SU{5)gut- In holomorphic 
gauge we choose 


{^xy) = AQi + d{X k)Q 2 + m{E^ + mxE^ ) myE2 ), (4.3) 

where Qi,Ef are E^ generators whose dehnition and all other details involving the Ug 
Lie algebra are given in appendix [A} Here A = fx^ibx — y), m,fh,yi and k are constants 
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with dimensions of mass and b, d are dimensionless constants. Notice that in terms of the 


background (3.12) we have chosen Ai = A and A 2 = —{d + 2)A — dn, with k being the 


distance between the two zeroes of these polynomials. As we will see, k will control the 
distance between the two Yukawa points pup and pdown of this model. 

As discussed in appendix [Bjthis background breaks SU (5)_l to S{U{2 )xU{2)xU{1))±, 
the representations of SU{5)± decomposing as 


SUi5)^ 

5 

10 


S{U{2) X U{2) X U{1))_ 


(4.4) 




2 

5 ’ 5 




(1, 1)6 _4 

^ ^ 5 ’ 5 


( 1 . 1 )- 


4 6 
5’5 


1,1 -2 

^ ^ 5’ 5 


( 2 , 1)1 _4 
^ ^ 5 ’ 5 


(1,2)_4 1 

^ ^ 5’5 


(2,2)1 1 

^ ^ 5’5 


where the subscripts denote the charges under the traces of the two U (2) factors. Each of 
these subsectors corresponds to a different 5 or 10 matter curve, and some of them will 
be associated to the 5'17(5)gut fields {10m, 5m, 5^/, 

In order to describe the different assignments of fields to matter curves and represen¬ 
tations of S{U{2) X U{2) X f/(l))± it proves useful to introduce a basis of the fundamental 
representation of 5 U^, that we denote by ei,... , 65 . The action of the background Higgs 
field on this basis reads 


/ 


= 


V 


Am 0 
m^x A 0 

0 0 d{X + k) 

0 0 

0 0 0 


m?y 


0 

0 

rh 

d(A -|- k) 
0 


0 

0 

0 

0 


\ 


(4.5) 


-2(1 -|- d)X — 2dK y 


in agreement with (3.13). A similar matrix $10 can be built for the antisymmetric rep¬ 


resentation of SU 5 , with elements of the form e* A Cj with i ^ j. From $5 and $10 and 


using eq.(3.2) one can detect the 10 and 5 matter curves, respectively. For this model 


one then finds three 10-curves and five 5-curves, each of them corresponding to a factor 
in the above decomposition of the 5 and 10 of SU{5)±, respectively. 

Assigning matter fields to matter curves as in the previous section is then equivalent 
to specifying their G_l = S{U{2) x 17(2) x 17(1 ))_l quantum numbers. One finds that for 


22 







the model with a satisfactory hierarchy of Yukawa couplings the assignment is 


10 


M ■ 



~ (2, 1)3 _2, 

/ ei A 65 

5m: 

5’ 5 

y 62 A 65 

~(1,1)_|^4, 

5d : (^3 A 64) 


(2,1). _i, 


(4.6) 


where we have also expressed these quantum numbers in terms of the basis ei,..., 65 , see 
appendix for further details. 

To cross-check this assignment let us read off each matter curve from the action of 
($) on the Gx quantum numbers of each matter field. More precisely, we define 
as the action of {^xy) on the part of (T^gut,”^^) C 248. We have that 

-(1 -|- 2d)\ — 2dK m 

rri^x —{l + 2d)\ — 2dK 

= —2A, <h|5^ = 2d{\ + k). 


$1 


10m 


1 ^ 

m \ ( 


5 *1* 5m = 

\ m X 

A / \ 


(4.7) 


The matter curves are then given by the vanishing of det <h| 7 e.Qy.j,, namely 


-m^x = 0, : {{1 + 2d)X + 2dK)f - m^x = 0 

^5c/ : A = 0, : d{X -|- k) = 0, 

in agreement with the discussion of the previous section and that of appendix 

The two types of Yukawa couplings are generated when these curves meet. In partic¬ 
ular, the up-type Yukawa 10 m x 10 m x 5u is developed at the intersection between 
and Ssy. On the other hand, the down-type coupling 10 m x 10m x 5u appears at the 
point where and Sg^ coincidej^ These points are 

Yu : n Sgy = {x = 1 / = 0} = Pup (4.9) 

^DjL • ^IOm ^ ^5m ^ ^0 ) y Vo} Pdown) 


with 




K 


Xq — n) 2/0 — ^ 1 + 






KbjY 

nY 


(4.10) 


We see that for each type of Yukawa there is a single intersection point and that, in general. 


these are not the same. It is the parameter k that controls the separation between them. 

^We are demanding that three curves in a surface meet, which does not look generic. However, due 
to gauge invariance, one of these equations is a linear combination of the others and the coupling is in 
fact generic. 
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For K = 0 both couplings are developed at the origin, whereas the two Yukawa points 
separate as |k| increases. 


Real gauge 

The previous background helds are in holomorphic gauge and we would now like to hnd 
a physical solution, namely one that satisfies the D-term equation. As explained earlier, 
this is achieved by performing an arbitrary complex gauge transformation and imposing 


(4.2), which translates into a set of differential equations for such transformation. 


More explicitly, consider the following transformation 


^ Ao^i + ig d g , 


(4.11) 


where g is an element of SU{5)^. We propose the following Ansatz 


g ^ gi(/Pi+/>2) 


(4.12) 


where Pk = [E^,Ej^] and /, / are real functions of {x,y). After the transformation, the 
Higgs and gauge fields are 

<h = XQi + d{X + k)Q 2 + m{e^Ef + mxe~^ E^) + fh{e^ E^ + mye~^ E^) (4.13a) 


Ao,,= -^{dfP, + dfP2). 


(4.13b) 


We can now plug these fields into the D-term (4.2) which yields equations for / and /. 
Taking the Kahler form as 

(4.14) 


oj = -{dx A dx + dy A dy), 


these become 


{d^ds: + dydy)f -w?\x\^e 

+ dydy)f - fh^\y\^e~‘^^). 


(4.15a) 

(4.15b) 


Following 22,26 , we take / to depend only on r^., the radial coordinate in the {x^x) 


plane. Then, eq.(4.15a) reduces to 


{ d?‘ IdX 1 


(4.16) 
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where s = ^{mrx)^, and the function h is dehned as 


= mrxe^’^ 


(4.17) 


This is a particular case of the Painleve III equation whose solution over the whole complex 


plane C can be found in 55 . However, since we are working in a patch of S'gut that 


contains the origin, we just need an approximate solution]^ which is 26 


f{rx) = logc + c + m rl ( 77 “ |^ ) + • • • 


(4.18) 


and similarly for /, replacing x —)■ y. In the previous equation the constant c needs to be 
hxed to the values 


= 3^/3 


£II1 

r[|] 


0.73, 


(4.19) 


if we ask for a regular solution for all values of r^. However, as mentioned above, we will 
not restrict to this particular choice since the actual value will be hxed only when all the 
global details of the background are specihed. 

4.2 Primitive fluxes 


The helds (4.13) dehne a consistent background that solves both the F and D-term equa¬ 


tions. We can still hnd a more general background by turning on additional gauge huxes, 
however, these cannot be generic since that would require modifying $. The most general 


hux that respects the Higgs held in (4.13) has to commute with it and be primitive on 


S'gut- If it also keeps the gauge group SU{5)gut unbroken such hux is of the form 
Fq = i{dx Adx — dy A dy){MiQi -\- M 2 Q 2 ) + f(dx Ady + dy A dx){NiQi + N 2 Q 2 ) ■ (4.20) 
As usual, the presence of such worldvolume hux is necessary to induce 4d chirality in the 


matter curves. The modes of opposite chirality 5, 5 and 10, 10 feel the background (4.13) 


in a similar way, and so whenever there is a zero mode solution for one chirality there will 
also be a solution for the opposite chirality. This is no longer true for the background 


hux (4.20), that will locally select modes of one chirality or the other depending on the 


®Far away from the origin the equations themselves receive corrections so it does not make sense to 
insist on solving the Painleve equation for the whole of C. 
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signs of Mi, Ni, i = 1,2. This chirality selection can be characterised in terms of a local 


chirality index 27 , as discussed in more detail in appendix [Cj 

Finally, following the standard strategy in F-theory GUTs, we break SU{5)gut down 
to the SM gauge group by turning on a flux along the hypercharge generator. Keeping 
the associated gauge boson massless amounts to imposing a global condition, which is 
invisible at the local level of our discussion. We then parametrise such flux as 


Fy = i 


Ny^dy Ady — dx A dx) + Ny{dx Ady + dy A dx) 


Qy 1 


(4.21) 


where the hypercharge generator is dehned as follows 

Qy = \ (^1 + 2^2 + 3 ^ 3 ) + \h^ . 
The total primitive flux is then 


(4.22) 


Fp = iQni^dy A dy — dx A dx) + iQs{dx Ady + dy A dx) 


(4.23) 


with 


Qr — NyQy — MiQi — M2Q21 Qs — NyQy + NiQi + N2Q2- 


(4.24) 


These fluxes will enter into the Dirac equation for the zero modes of our model. As a result, 
each of the MSSM chiral zero modes will feel a different flux depending on their quantum 
numbers (and in particular its hypercharge) and will then develop a different wavefunction 
prohle. As mentioned before, these flux differences will not affect the Yukawas at the 
holomorphic level, but they will enter into the hnal expression for the physical Yukawas. 
We have then gathered the flux felt by the different MSSM sectors of the present model 
in table [TJ and in particular the effective combination of fluxes qr and qs that will be 
crucial for the computations of section 


4.3 Residue formula for Yukawa couplings 

The computation of the holomorphic Yukawa couplings can be performed via dimensional 
reduction of the 7-brane superpotential 

W = ml [ Tr (<F A F) + -^oTr {F A F) . (4.25) 

Js 2 
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MSSM 

Sector 

S{U{2) X U{2) X U{1))± 

Gmssm 

Qr 

qs 

Q 

10m 

(2, 1)3 _2 

(3.2)_. 

-liVy-Mi 

MNy + Ni 

u 

10m 

(2, 1)3 _2 

5 ’ 5 

(3,1)| 

INy - Ml 

INy + Ni 

E 

10m 

(2, 1)3 _2 

5 ’ 5 

(l.l)-i 

-Ny - Ml 

-Ny + Ni 

D 

5m 

( 2 , 1)1 _4 

5 ’ 5 

(3.1). 

— |iVy -|- Ml 2 M 2 

-\Ny + Ni-2N2 

L 

5m 

( 2 , 1)1 _4 

(1.2). 

|iVy -|- Ml 2 M 2 

\Ny + Ni- 2 N 2 

Hu 

5c/ 

(1,1)_6 4 

(1.2)_. 

-liVy + 2Mi 

-\Ny-2Ni 

Hd 

5d 

(1,1)_4 6 

^ ^ 5’5 

(1.2). 

liVy - 2 M 2 

\Ny + 2 N 2 


Table 1: Different sectors and charges for the model of this section. Here and qs 


are the operators (4.24) evaluated at each different sector. All the multiplets in the 


table have the same chirality. 


As discussed in section the second term in (4.25) is due to the presence of non- 
perturbative effects in the compactihcation. In particular we have that 9q is a holomorphic 
section on S and e is a parameter that measures the strength of the non-perturbative ef¬ 
fect. We note that the presence of this additional term will eventually change the BPS 
equations for the background that we previously solved, and so the background prohles 
for $ and E will have 0(e) corrections. One can however show that these corrections do 


not affect the computation of holomorphic Yukawas 25 , and so they can be ignored in 
what follows 13 

The zero mode equations can be derived by expanding the 7-brane fields into back¬ 
ground and fluctuations 

<h = ($)-1-, A={A) + a, (4.26) 


and linearising in fluctuations the F-term equations derived from (4.25). We obtain 


d{A)a = 0 , 

d(^A)V> = i[a, (<h)] - edOo A (9(^)0 -f 9(^)0^) • 


(4.27) 


^The full superpotential expression involves terms of the form 9kSTi{^^yF^), k > 2, but suppressed 


by higher powers in m* 20 24 25 . Hence their contributions will be less relevant than those from 00- 
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A similar procedure can be applied to the D-term equation, but since we are simply 
looking at the holomorphic part of the Yukawa couplings in this section we will postpone 
that discussion to the following section. It is possible to solve explicitly for the system 


(4.27), the solution being 


(4.28) 


a = d{A)^, 

ip = h- i[(<h), + ed9o A (a^ - d{A)0 , 

where is a section of ® ad(ii^8) and h is a holomorphic section of 


ad(i?8)- We stress that while the solution (4.28) contains dependence on a'^ which may 


in principle introduce some non-holomorphic terms in the 4d superpotential, these terms 
will appear only in total derivatives and so they will not eventually appear in the resulting 


4d superpotential 25 . Using this solution it is possible to prove 17,22,25,26 that the 


Yukawa couplings are 


.mZ 


Y = -i^ / Tr(hAa(^)eAa(^)0. 


(4.29) 


It is also possible to write the Yukawa couplings as a residue evaluated at the Yukawa 


point. We simply quote here the result referring to 25,26 for the general proof: 


Y = m^TrVabc Resp [r]°‘r]’'h^y\ = mln'^fabc j il^rj'^Kydx A dy , 


(4.30) 


where C can be continuously contracted to a product of unit circles surrounding the 
Yukawa point p without encountering singularities in the integrand and we dehned the 
function p as 


p ^hxy 'i^d)x0{)dy hxy^ ladyBoOx 


(4.31) 


4.4 Holomorphic Yukawa couplings for the Eg model 

Let us hnally discuss the structure of the Yukawa couplings that arise in the present local 
Eg model. We focus our attention on the sector involving only the MSSM helds so we will 
have only two Yukawa matrices, namely the x 10 m x and the 10m x 5m x 5h- 
The functions hxy for the different helds are 

hioM = Yo,i ml-\bx - yf-^ = 75 ,^ ml-\h{x - xq) - {y - yo)f~^ 

h^H = lu h^H = 7 d , 
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where (xq, yo) corresponds to the coordinates (4.10) of the down-type Yukawa point Pdown, 
the constants 710,i, 75,i, 7(7, Id are normalisation factors to be hxed in the next section and 


i = 1, 2, 3 is a family index. Using these we can compute the functions rj in (4.31) which 


in turn are needed to compute the holomorphic couplings via the residue formula (4.30). 
Such 7 ’s are computed in appendix where the following Yukawa couplings are found: 
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7Z. 710,175,2 
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^ 2dp^75,3 
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d J 2p^7io,3 
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(4.33) 


e — + bOy), 


K 


K = 


m* 


m 

Pm — 7 ) 

mi 


Pp 


Pm710,3 

1 

mi 


(4.34) 


and where we have kept terms linear in e. Notice that Y[/ is of the form (3.14) and Yd/^. 


reduces to (3.16) in the limit k —)■ 0 . A non-vanishing k distorts the form of but it 

does not spoil its hierarchical structure of eigenvalues. In fact, as we will see in section]^ 
K will only enter in the CKM matrix describing quark mixing angles. 


5 Zero mode wavefunctions 

So far our discussion of the Yukawa couplings has been restricted to the holomorphic level, 
namely we have been discussing the cubic couplings that appear in the four dimensional 
superpotential. The missing ingredient in the computation of the physical couplings 
involves the normalisation of the wavefunctions as well as their kinetic mixing. Indeed, 
recall that the Yukawa couplings that are measured experimentally can be compared 
to those that appear in the Lagrangian once the kinetic terms for the chiral helds are 
canonical. In order to compute such kinetic terms, we necessarily need to solve the 
equations of motion for the zero modes in a real gauge, which introduces the dependence 
on the worldvolume flux densities. 

In this section we will solve for such zero mode wavefunctions of the Ug model of 
section and compute the kinetic terms for them. Solving analytically for real gauge 
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wavefunctions is a much more complicated problem than doing it at the holomorphic level, 
especially for T-brane models like ours, and the problem becomes particularly involved 


when non-perturbative corrections are taken into account. However, similarly to 26 these 


zero mode equations can be solved for a certain region of parameters of the local model, 
allowing to see how such wavefunctions depend on flux densities. Our approach will be 
to hrst consider the perturbative case and compute the kinetic terms, where no kinetic 
mixing arises for the choice of wavefunction for each family that we make. Second, we 
include the non-perturbative corrections and argue that they do not change the result. 
Many computational details will be relegated to appendix]^ (see also 24-^). 

That is, in this section we will see that i) there is no kinetic mixing between families 
and ii) non-perturbative corrections do not affect their kinetic terms, at least at the level 
of approximation that we are working. As a result, to compute physical Yukawas one may 
combine the residue computation of the non-perturbative holomorphic couplings with the 
computation of the normalisation factors 710 , 75 at tree-level that we do in the following. 


and which should be inserted in eq.(4.33) to obtain the physical Yukawas 


10 


5.1 Wavefunctions in the perturbative limit 

Before including non-perturbative corrections, the equations for the zero modes are ob¬ 


tained from (4.1) and (4.2) by expanding to linear order around a given background as in 


(4.26) which yields 


d{A)a = 0 , 
di^A)V (^)], 
wAa(A)a = ^[($),(p]. 


(5.1a) 

(5.1b) 

(5.1c) 


where (<h) and (A) correspond to the background in real gauge. These can be solved for 


every particular sector using the techniques in 24-26 and in the following we just quote 


the result for the relevant sectors, namely those that appear in table The details of the 
computation can be found in appendix 

Alternatively, one may directly compute the physical Yukawas by performing the triple overlap of 


real gauge zero modes. Similarly to 25 26 one can show that both approaches give the same result. 
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We use the following notation for the zero modes, 





al 




p,s 


Ky y 


(5.2) 


where Ep^g denotes the particular set of roots, labeled by s, for each sector p. In the case 
of the up and down-type Higgses, s only takes one value since these are only charged 
under an Abelian subgroup of S{U{2) x U{2) x U{1)). On the other hand, the matter 
sectors transform as doublets of the hrst U{2) factor so s takes two different values in 
that case. 


Higgs wavefunctions 


The solution for the 5u sector is 


= 


/ —jAl \ 


7c/ 


■ (u-\u 


V 


Xu E, 


u 


(5.3) 


/ 


where 




Xuix,y) = e 2 

and Xu is a function of the flux densities and intersection parameters given as the lowest 


solution to the cubic equation (E.6). Also, Cu = ^2(Xu-^s) 


_ ^ui^u-qR-gs) 


Similarly, the solution for the 5^^ is 


— 


/ \ 


7d 


■ C,d-\d 


e - xo,y - yo) Ed 


(5.5) 




with 


XD{x,y) = 62 


— e^(.\h^-\y\^)-qs(^y+y^)+ix-y)(CD^-{^D-CD)y) 


(6.6) 


and where ijj is defined in (E.IO). Finally, Xd is the lowest solution to (E.13) and (d = 

>^D{>^D-qR-qs) 

‘^{^D-qs) 
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Matter wavefunctions 


These sectors are a bit more involved because the helds are charged under the T-brane 
background. Given our choice of background, both the 10 m and 5 m transform as doublets 
of the hrst SU{2) factor in the decomposition of SU{5)^. Thus, we write the solution as 




Et + 


\ K / 


El — _E I , 


(5.7) 


\ ) 

where we use a + to denote the upper component of the U{2)i doublet and — to denote 
the lower one. The zero mode equations in these two sectors turn out to be rather 
complicated to solve in general, but it is still possible to hnd approximate solutions in the 
limit /i, K <C m. The real wavefunction for the 10 m is 

7 ^Aio \ /o\ 


^■io = 


7io 


AioCio 


V 


0 


/ 


0 

w 




(5.8) 


where Aio is the negative solution to the cubic (E.25) and Cio = ~?s/(-^io ~ Qr)- Finally 
the wavefunctions Xio are 


_ g^(hk-|yk)-9sG?/+J/2:)+Aioa:(a;- 


+ Cloa:), 


(5.9) 


where are holomorphic functions oi y + (Ciox and i = 1,2, 3 is a generation index. As 


m 


25,26 we choose these holomorphic functions in the following way 

Aaiy + Cioi) = ml-'iy + Cioa:)““‘. 


(5.10) 


The solution to the 5 m is very similar and reads 

( 1^5 \ / 0 ^ 

e'*+f/\^{x,y-,,/a)E++yi 




-i^ 


\ 


0 


/ 


0 

VV 


f/‘^j^^(^x,y-v/a)E^ (5.11) 


with 'tp dehned in (E.29). Also, As is a function of the fluxes and intersection parameters 


dehned as the lowest solution to (E.31) and Cs = —^^/(As —gj?). Finally the wavefunctions 
Xs are 

y) = ^ (512) 
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where g\ are holomorphic functions of y + C,^x and i = 1,2,3 is a generation index. 
Analogously, the family functions are 

9l{y + = ml-\y + • (5.13) 


5.2 Normalisation factors 


Once we have the perturbative wavefunctions we can compute the normalisation factors 
and kinetic mixing. The appropriate scalar product is given by 


K'’ = = mt / iv(y‘'.yj)dvoi, 

Js 

as can be seen by performing the dimensional reduction. 


(5.14) 


Given the choice of family functions (5.10) and (5.13) we hnd that the kinetic terms 


(5.14) are diagonal, so we only need to compute the corresponding normalisation factors. 
We hnd 


i7^r 

i7io.r 

l75.f 


4 

^2 


(2C[/ + + 2Cf/ — 2A^) + {qs + Xu)"^ 

Ti* / 4/i^ + + iCu — 

h {‘^Cd + Qr){9r + ~ 2A_d) + {qs + XpY 

) 4^2^4 + (2^ + [C,D - Ad)^ 

c 1 


ml7r‘^{3 — j)\ _1_L ^^^10 _1_ 

2Aio+iJfl(l+CiQ)-»n2c2 m4 2Aio+(jfl(l+C^g)+m2c2 



C 1 

m%2(3 — j)!_1_L ^^^5 _1_ 

2A5+<Ji?{l+C5)-"i2c2 m* 2\5+qji{l+C^)+m'^c^ 



(5.15a) 

(5.15b) 

(5.15c) 

(5.15d) 


Recall that the parameters Xp and (p for a given sector p, depend on the hux densities felt 
by such a sector and, in particular, depend on the hypercharge hux. Thus, each MSSM 
multiplet in a given GUT multiplet will have diherent normalisation factors. 


5.3 Non-perturbative corrections to the wavefunctions 

The computation of the kinetic terms performed above allows to obtain the physical 
Yukawa couplings at tree level. However, since we are interested in the leading corrections 
induced by the non-perturbative ehects, we need to compute the normalisation factors 
and mixings at 0{e). To do so, we will solve for the zero mode equations in real gauge 
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including all 0(e) corrections. As we will see, it turns out that no mixing is generated and 
the normalisation factors are not corrected at this order, so we may use the ones obtained 
earlier. 


Following section 4.3, the F-term equations for the zero modes at 0{e) read 


d{A)a = 0 , 

d(^A)V> = i[a, ($)] - edOo A (9(^)0 + d(^A)a^) • 


(5.16) 


which have to be solved together with the D-term equation (5.1c) that remains unchanged 


1^ . As in the tree-level case, we quote the relevant results for each sector and relegate 
the computations to appendix We start with the Higgs sectors that are not charged 
under the T-brane and then consider the more involved case of the matter sectors. 


Higgs sectors 

The solution to the non-perturbative zero mode equations for the sector 5{/ is 


/ 


= 


lu 


4 (u 

: iCu-^u) 


\ 


^np _ -\y\'^)-qs(xy+yx)+{x-y){(4ux-(\u-Cu)y)) 


\ 1 y 

The 0{e) non-perturbative correction is 


(5.17) 


<^2 


= -^{Cux-iXu-Cu)yY{0x; + 9y) + ^{x-yf + ^{x-y){Cuy + {Xu-Cu)x) (5.18) 


with the constants 5i, 62 given by (E.42) and (E.43) respectively. The solution to the 5^ 


is essentially the same and can be obtained by performing the replacements explained in 
appendix so we do not write it explicitly. 

Now one can see that this particular correction to the wavefunction will not generate 
a correction to the normalisation factor at order e. The reason is that the extra terms 
that appear in the integrand of (5.14) will be those in ( |5.18[ ) and its complex conjugate 
which are not invariant under the rotation {x,y) —)■ e*“(a:, ?/). 
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Matter sector 


As shown in the appendix, the strnctnre of the solntion for the 10 m sector is 

/o\ /o^ 

f O(e^) = 


"^10+ ~ 


+ e 


vv 


0 

V/ 


0 

V/ 


+ 0(e^). (5.19) 




and similarly for the 5 m- This strnctnre already shows that the 0(e) corrections to 
the kinetic terms of the matter sectors vanish, even withont specifying their explicit form. 


Indeed, from (5.14) snch corrections will be proportional to the scalar prodncts 


and • "^[ 0 +’ "'^here the superscript (0) denotes the tree-level term and (1) the 0 ( 


10- 

e) 


correction. Given the solution (5.19), we see that those products are trivially zero. 


6 Physical Yukawas and hierarchies 

Combining the results of the last two sections one hnds the following physical Yukawas 
for quarks and charged leptons in our local Bg model 
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with the dimensionless complex parameters dehned by 
e = e(dx + K = 


K 


m* 


m 

Pm — 7: 

mi 


- 

mi 


(6.1a) 


0(6^ 


(6.1b) 


0 {t 


(6.1c) 


( 6 . 2 ) 
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and the normalisation factors 'yfoj, 75 ^ given by (5.15). The snperscript a denotes the 
particnlar MSSM chiral mnltiplet within 10 m or 5 m, namely the hrst colnmn in table 
We wonld like to see if this strnctnre can reprodnce empirical data for charged fermion 
masses. Since onr resnlts apply at the GUT scale, presnmably at aronnd 10^® GeV, ex¬ 
perimental valnes at weak scale need to be rnn nsing the renormalisation gronp eqnations. 
Table shows the extrapolation to the nnification scale of snch observed qnantities taken 
from in the context of the MSSM. These depend on the parameter tan/? that con¬ 
trols the relative magnitnde of the vevs of Hu and Ho- In particnlar, we have that 
mr,b = y^bV cos/3 and rrit = YtV sin/3 with V = \JV^ -f Vj 174 GeV. In the following 
we will discnss the comparison between the experimental resnlts with the predictions from 
onr local & model. 


tan/3 

10 

38 

50 

rriu/mc 

2.7 ±0.6 X 10-3 

2.7 ±0.6 X 10-3 

2.7 ±0.6 X 10-3 

mjmt 

2.5 ±0.2 X 10-3 

2.4 ±0.2 X 10-3 

2.3 ±0.2 X 10-3 

md/rris 

5.1 ±0.7 X 10-2 

5.1 ±0.7 X 10-2 

5.1 ± 0.7 X 10-2 

rris/mb 

1.9 ±0.2 X 10-2 

1.7 ±0.2 X 10-2 

1.6 ±0.2 X 10-2 

rrie/m^ 

4.8 ±0.2 X 10-3 

4.8 ±0.2 X 10-3 

4.8 ±0.2 X 10-3 

m^/m-r 

5.9 ±0.2 X 10-2 

5.4 ±0.2 X 10-2 

5.0 ±0.2 X 10-2 

Yr 

0.070 ±0.003 

0.32 ±0.02 

0.51 ±0.04 

Yb 

0.051 ±0.002 

0.23 ±0.01 

0.37 ±0.02 

Yt 

0.48 ±0.02 

0.49 ±0.02 

0.51 ±0.04 


Table 2: Running mass ratios of leptons and quarks at the unification scale from ref. 56 
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6.1 Fermion masses 


The masses for quarks and charged leptons will directly depend on the eigenvalues of the 


physical Yukawa matrices. From (6.1) we see that such eigenvalues are 


Yt 

Yb 

Yr 


Y lU 7?o,37io,3 

2pmP/i 
2 Q D 

^ 7k),37^3 

2,dpmpfi 

-yp T'fo,371:3 

‘Idpmpfj, 


Y — 2 'yu •T'?0,27io,2 

4pmp2 ) 

2 Q D 

Y T'iO,2T'^2 

^ 4d2p^p2 , 

Y — - 70 7fo,275,2 


= O(e') 

Id = O(e') 

n = 0(6-2), 


(6.3) 


which makes manifest the mass hierarchy (0(1), 0(e), 0(6^)) between families. However, 
it still remains to see if the data of table can be reproduced via these expressions, 
and if affirmative for which range of values for tan 13. This question is non-trivial in the 
sense that the normalisation factors 710 , 75 are complicated functions of the multiple flux 
densities present in the model, which makes it hard to proceed analytically]^ 


Mass ratios 


While the Yukawa eigenvalues are complicated functions of the local flux densities, fermion 

In 


mass ratios have a much simpler dependence on them, as already noticed in 25 


particular, let us consider those mass ratios between the second and third generation that 
are independent of tan (3. These are 


rrtc 
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‘3‘dpfi 
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El 

T 
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INy 

‘^dp^ 
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(6.4a) 


(6.4b) 


(6.4c) 


where g®, are the linear combinations of flux densities that appear in table and we 
have dehned the quotients 


Ml Ml + 2M2 

X = —^ y = -=- 

Ny Ny 

^^Notice that when embedded into a global model, this large number of flux densities should depend 
on a few Kahler moduli of the compactihcation, and in this sense many of the free parameters that are 
present in the local approach become constrained. 


(6.5) 
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We see that for these mass quotients the dependence on the normalisation factors for the 
Higgses drops and we obtain fairly simple formulae in terms of a few flux densities. More 
precisely, besides x and y these three ratios depend on two more parameters, namely 
\d\ and |eiVy/2p^|. Furthermore, by considering quotients of ratios we can eliminate the 
dependence on the last parameter, since 


mjmt 

ms/rrib 

m^/rrir 

rris/mb 


= \d\ 


= 3 


^+3 

y-l 


1/2 


(x - 1) (l/ + \) 


1/2 


(6.6a) 

(6.6b) 


(3x- I) (3?/- 1) 

So one may proceed to constrain these three parameters of the local model in terms of 
two empirical quantities, which from the data of table read 


rric/mt 


ms/mh 


exp. 


m^/mr 

nis/mi, 


exp. 


0.13 ±0.03 

3.3 ± 1. 


(6.7a) 

(6.7b) 


Finally, recall that the values of x, y are constrained from the results of appendix For 
Ny < 0 we have that x < —2/3, y < —1/2, while for Ny > 0 we have that x > 1, i/ > 1/3. 
We hnd that it is easier to £t the above empirical values for the latter case and by taking 
a small value for \d\, as illustrated in figures and For instance, taking d ~ 0.02 we 


can use (6.6) together with (6.7) to estimate x,y dehned above, namely 


X =5±3 
1/ =0.45 ±0.05. 


( 6 . 8 ) 


Then, using these approximate values for the fluxes we can hnd the order of magnitude 


of leNy/p^l by htting one of the mass ratios in (6.4), which yields (see hgurep 


eNy 

Pm 


(1.3 ±0.7) • 10"’. 


(6.9) 


Notice that, up to now, we did not specify the value of any hux density but only quotients 
of huxes. In the following we discuss the absolute value of Yukawa couplings for which 
we actually need the hux density values. 
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Figure 1: Ratio of ratios (6.6a) where the horizontal axis is x = —Mi/Ny and the vertical 
is y = (Ml + 2 M 2 )/Ny- The right figure displays the value for \d\ = 0.1 and the left one for 


\d\ = 0.02. In the latter, the red dot corresponds to the values chosen in (6.12). 


Yukawas for the third generation 


So far we have discussed ratios of Yukawa couplings for which we have simple expressions. 
However, for the Yukawa couplings themselves, we do not have such simple results and 
it is much harder to understand how these depend on the parameters of the model. In 
particular, these depend on flux densities and not just on quotients so let us estimate 
their order of magnitude. 

The fluxes that are not along the hypercharge generator, i.e. Mi, M 2 , Ni, N 2 , deter¬ 
mine the number of chiral multiplets for the matter sectors. On the other hand, the 
hypercharge fluxes Ny and Ny are also subject to quantisation conditions since they are 
responsible for the doublet-triplet splitting. Thus, all of the flux densities satisfy that 
I Jj, F| ~ 27r, where S is a matter curve and if we take F to be approximately constant we 
find that |F| ~ 2tt/V y, with Vy the volume of E. Furthermore, we know that the volume 
of the GUT divisor, Vgut, is related to the coupling constant acur as (see e.g. [^) 


acuT — 


27r^gf^ 
^tt^GUT ’ 


( 6 . 10 ) 


where the string scale rrist is related to m* by = {27i)^gsm‘^. From unihcation of the 
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Figure 3: Ratio ms/rrih that fixes the order of magnitude of |eiVy/2p^|. 

1 /2 

gauge couplings, we expect acuT — 1/24, so assuming Vs — ^gut^ 


|F| 

1 ^ 1 

f ‘^OLguT^ 

0.3 


\ 9s ) 

' ^ sT- 


( 6 . 11 ) 

which tells us the order of magnitude for the the flux densities. Notice that our approach 
relies on having diluted fluxes (as well as intersection slopes) so the coupling constant 
should not be arbitrarily small. 


Now one can perform a scan for flux densities in the ballpark of (6.11) and see whether 
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there is a region that allows to ht the Yukawa couplings for the third generation of U and 
D-type quarks as well as charged leptons. We hnd that this is indeed possible and, for 
instance, taking the following flux densities (in units of 

(Ml, M 2 , Ni, N 2 , Ny, Ny) = (-0.16, 0.09, -0.501, 0.501, 0.03, -1.0), (6.12) 

together with intersection angles 


(p™, p^, c, d) = (0.09,0.004,0.53, 0.025), 


(6.13) 


we hnd 

Yt = 0.46 W = 0.08 = 0.15, (6.14) 

consistent with the experimental values in tablefor tan/? ~ 10 — 20. Also, using these 


parameters we get from (6.9) that 


e = 1.3-10 


-4 


(6.15) 


which is small enough to be consistent with its non-perturbative nature. 


-0.4 


- 0.6 


Ny - 0-8 


-1.0 


-1.2 

0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 

Ny 



Figure 4: Ratio of r and b mass. The red dot corresponds to the values chosen in (6.12). 


As we can see, our model allows to accommodate a large top Yukawa coupling, which 
is usually troublesome in perturbative type II GUTsj^ Also, the hypercharge hux may 


“^Notice that the expression for the top Yukawa in our model is exactly like in the Eq model of 26 . 
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induce the correct difference between Yf, and YV- Recall that these come from the same 
Yukawa in the SU{5) GUT but, due to the hypercharge breaking, they are different even 


at Mgut- From (6.3) we have that 




7 ^, 375(3 

7^,375^3 


5 


(6.16) 


which, for vanishing hypercharge fluxes, is exactly one. However, for our particular choice 
of parameters we hnd that 

7 = 1.81, (6.17) 

lb 

consistent with the observed ratio. 


6.2 Quark mixing angles 

Let us now analyse the quark mixing angles for this model. Recall that the CKM matrix 
is dehned in terms of the unitary matrices Vu and Vo such that they diagonalise the 
Hermitian product of quark Yukawa matrices YYb More precisely we have that 

Mu =VuYuY{jV{j (6.18a) 

Md=VdYdYIvI (6.18b) 


with Mu and Mu diagonal. We then define the CKM matrix as 


VcKM — VuV^. 


(6.19) 


Directly applying these dehnitions to (6.1a) and (6.1b) and not taking into account 
their 0 {e‘^) corrections, one hnds that the following matrices satisfy (6.18) 

/ 1 n \ 


Vu = 


1 

0 

* Q 
7 io,i 


\ 2p*7]^Q 3 

/ 


0 

1 

0 


- Q 


0 

1 


Ud = 


! 710,1710,2 

2 kY\ 7 ( 0,1 
Pm 


'iPMlPMp(7(o,3)^ 


2PM7fb,3 


2|PmP(7?o,3)" 

^* 7 ( 0,2 

Pm 7 io ,3 


2 kl 

Pt^ 

~ Q 
^7io,2 

Pm7?o,3 


Q 

7ioi 

2Pm7(o,3 




|feP(7l^0,2) 

2|Pm 177 ( 0 , 3 )^ / 


(6.20a) 


(6.20b) 
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Taking into account O(e^) corrections would modify these expressions, in particular those 
related to the rotation angles for the hrst family. In particular, one expects that the hnal 
rotation matrices are to a good approximation of the form 


Vcj = RuVu, 


Vn = RnV, 


DVD 


( 6 . 21 ) 


with 


/ 


Ru,D ^ 


OiU,D e 

.2 1 


2 0 
0 




V 


( 6 . 22 ) 


/ 


—au,D e 

0 0 1 

and where au, an are 0{1) unknown rotation angles. These extra rotations will modify 
the value of several CKM matrices elements, but leave untouched the mixing between the 
top and bottom quarks. To the degree of approximation that we are working, we hnd 
that such entry reads 


I«P(7?0,2) 


\Vtb\^l- 


2|p.P(7?o,3)^ 

Since the experimental value for this quantity is 


= 1 - 


^ gi?,Q 

^Pl ' 


(6.23) 


|htb|exp. — 0.9991 


(6.24) 


we hnd a typical value |«)| ~ 10 ^ — 10 In particular using the values (6.12) we hnd 


\k\ = 2.7-10 


-3 


(6.25) 


That justihes the approximation k m made in appendix 1 Interestingly, this result 
has a direct geometrical interpretation. We have that the mixing between the second and 
third family is roughly given by 


\/l - \Vtb\ ^ 







oc m^\bxo - yo\ 


(6.26) 


where (xq, yo) are the coordinates of separation of the Yukawa point Pdown with respect 


to Pup, see (4.10). Hence we hnd that the mixing between the second and third family is 
proportional to the separation of the two Yukawa points along the particular direction bx— 
y. This combination of x and y is nothing but the complex variable that the holomorphic 


piece of the matter wavefunctions depend on, c.f. (4.32). Hence, this mixing ehect can 
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be understood as a change of wavefunction basis when moving from one Yukawa point to 


the other, in agreement with the results of 15 . Notice that the experimental value (6.24) 
then translates into a separation between Yukawa points which is roughly so 


Pup and pdown need to be relatively close to each other, as anticipated in 13 


The other two mixing angles are more difficult to estimate, as this would involve ob¬ 
taining explicit expressions for Ru and Rd. In any case, it is unlikely that a rotation 
of this order in e will generate the large experimental value for the Cabibbo angle. At 
this point one should however recall that in our analysis we have made certain approxi¬ 
mations related to the fact that we are describing our system in a small neighbourhood 


of the Yukawa points. These approximations include taking a flat metric (c.f. (4.14)), 
the limit /i m that neglects the curvature of the matter curves, taking Ai and A 2 in 


(3.12) as simple functions, and taking the holomorphic piece of the matter wavefunctions 


as monomials (c.f. (4.32)). These approximations are justihed for analysing the physical 
couplings that involve the two heaviest families, as their real wavefunctions are by con¬ 
struction localised near the Yukawa points. But this need not be so for the hrst family, 
which in specihc constructions may not even be described as a local chiral mode in the 


sense of 27 . Hence, we expect that the mixing angles involving the hrst family will be 
particularly sensitive to curvature corrections of the matter curves and S'gut- Finally, 
following the arguments in [^, one may estimate that corrections effects of the order 
precisely those necessary to generate a realistic Cabibbo angle. 


7 Conclusions 

In this paper we have considered the structure of both up and down-type Yukawa couplings 
in local F-theory SU{5) GUTs that are generated at a region of Eg enhancement. In 
particular, we have analysed several ultra-local models that realise the breaking to SU (5) 
via reconstructible T-branes and have rank one Yukawa matrices. We have then included 
the effect of non-perturbative dynamics (Euclidean D3-branes or gaugino condensates 
on 7-branes) and studied its impact on the Favour structure of our local models. In 
particular, we have seen that all families of quarks and charged leptons become massive 
once that these non-perturbative corrections are taken into account. 
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The hierarchical structure that results from combining tree-level and non-perturbative 
effects is already manifest at the level of the holomorphic Yukawas, in terms of the strength 
of the non-perturbative effect e which we use as an expansion parameter in onr computa¬ 
tions. This allows to classify the whole set of models under stndy by the eigenvalue hier¬ 
archy present in the np and down-type holomorphic Yukawa couplings. We have identified 
as promising models those whose eigenvalne hierarchy is of the form (0(1), 0(e), 0(e^)) 
for both types of holomorphic Yukawas, as this structure has proven to be successfnl in 


reprodncing empirical fermion masses for reasonable valnes of e in previous work 25,26 


We have carried such classihcation of models in section dnbbing each of the models 
in terms of the block-diagonal structure of the T-brane profile for the 7-brane held <h. We 
have seen that the 4-|-l model has the appropriate hierarchy of Yukawa couplings, bnt has 
the less attractive featnre of having both Higgses in the same matter curve so one would 
quite likely hnd a large /i-term in a global completion. This problem can be easily solved 
in any of the fonr 3+2 models analysed snbsequently, as they contain several 5-matter 
curves. However, we hnd that two of these models have vanishing up-type Yukawas at tree 


level, in a similar fashion to the vanishing resnlt of 23 . It would be interesting to acqnire 
a better nnderstanding of these vanishing results. The two remaining 3-|-2 models have 
rank-one tree-level Yukawas but fail to generate an appropriate hierarchy when inclnding 
the non-pertnrbative ehects. Finally, we have analysed fonr diherent 2-|-2-|-l and found 
that two of them show the desired pattern of fermion masses, and one in particular exhibits 
interesting mechanisms to generate realistic nentrino masses and /i-term. 

Consequently, we have performed a more detailed analysis this last 2-|-2-|-l model in 
order to compute its physical Yukawas. We have specihed the whole set of worldvolume 
fluxes that account for chirality, SU{5) breaking and doublet-triplet splitting. We have 
also included the non-primitive fluxes that solve for the T-brane backgronnd eqnations 
of motion. Finally, we have compnted the wavefnnctions for the matter and Higgs sec¬ 
tors in real gange at leading order in the non-perturbative parameter and extracted the 
normalisation factors for the 4d matter fields. 

In section we have combined the holomorphic Yukawas with the matter field nor¬ 
malisation factors, obtaining the physical Yukawa couplings. While the former does not 
depend on the fluxes, the latter does introduce a dependence on the fluxes which has 
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drastic phenomenological implications. Indeed, as we have seen these are crucial to £t the 
observed fermion masses and in particular to generate a difference between the D-type 
quarks and leptons mass ratios. Typical values of flux densities allow to obtain a large 
top quark Yukawa, with e 10 ^ and tan {3 10 — 20. Finally, we have considered the 

possibility in which both Yukawas are realised at different points within the region of 
enhancement. We have seen how this separation translates into the CKM matrix and 
checked explicitly that the distance between the points generates quark kinetic mixing. 
Thus, as expected, the fact that the CKM is mostly diagonal in encoded in the fact that 
the two Yukawa points are very close to each other. 

The above results are very promising for the F-theory GUT programme, and it would 
be interesting to extend the present analysis in several directions. One obvious direction is 
to perform our analysis beyond the leading order in perturbation theory in e. This would 
allow to gain some control over the lightest family of chiral matter, and in particular 
to analyse in more detail the CKM matrix. Nevertheless, as argued in the last section 
full control over the lightest family cannot be achieved until we extend the description 
of our model to a region containing the matter curves that host the SU(5) multiplets 
and IQm) and include curvature corrections to our computations of the wavefunction 
normalisation factors. In general, it would be extremely interesting to promote our ultra¬ 
local model to a local one where S'gut is a compact four-cycle. This would allow to 
understand the plethora of flux densities that enter as free parameters in our model 
in terms of a few Kahler moduli, and see if the flux relations that we have obtained 
are compatible with the geometry of Ygut- From a broader perspective, it would be 
important to embed our model within a class of fully-fledged F-theory compactihcations, 
and interpret each of the local parameters in terms of complex and Kahler moduli of 
the compactification. One may then see whether the region of parameter values that are 
needed to reproduce the flavour structure observed experimentally can indeed be reached 
when scanning through the F-theory landscape. 
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A machinery 

The Lie algebra of Eg consists of 248 generators Q^- We will work in the Gartan-Weyl basis 
{Hi, Ep} of Cg and where the generators with i = 1,..., 8 form a basis of the Gartan 
subalgebra and the remaining 240 roots are chosen to satisfy the following commutation 
relations 

[Hi, Ep] = PiEp . (A.l) 

This allows to represent the roots with a vector of charges under the Gartan subalgebra 
and for the case of eg the roots are 

( ±1,±1, 0,0,0, 0,0,0 ), (^±i, ±i, ±i, ±i, ±i, ±^, ±^, ±0 with even + . (A.2) 

For our purposes we need to decompose the Eg Lie algebra as Eg —)■ SU {5)gut x SU (5)_l. 
In particular the branching rule for the adjoint representation of Eg is the following 

248 ^ (24, 1) © (1, 24) © ((10, 5) © c.c.) © ((5,10) © c.c.). (A.3) 

We identify the roots in the adjoint representation of SU{5)gut as 

(0,0,0,1,-1,0, 0,0), (A.4) 
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which together with the Cartan generators: 


= H2 = H5-He, H3 = He-Hr, = Hr - Hs. (A.5) 

give the adjoint representation of SU{5)gut- The adjoint of SU{5)± is consists of the 
following roots 


(± 1 , ± 1 , 0 , 0 , 0 , 0 , 0 , 0 ) ± ( 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 



1 

1 

1 

1 

1 

1 

1 , / 

i 2’ 

“2’' 

“2’' 

“2’' 

“2’' 

“2’' 

■2'“2j + l 


1 1 111111 
2 ’ 2 ’~ 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 
11111111 
2’2’~2’~2’~2’~2’~2’~2 


(A. 6 ) 


and Cartan generators 




(A.7) 


where H^l = We will label the roots of the adjoint of SU{5)± as follows 

ii;± = ±(0,-1,1,0,0,0,0,0), 


E^ = ± 


11111111 

' 2 ’ 2 ’ 2 ’~ 2 ’~ 2 ’~ 2 ’~ 2’~2 


^ 3 ^ = ±(- 1 , 0 , 1 , 0 , 0 , 0 , 0 , 0 ), 

= ±(- 1 ,- 1 , 0 , 0 , 0 , 0 , 0 , 0 ), 


Ei = ± 


11 111111 
2 ’ 2 ’ ~ 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 


^ 6 ^ = ±(- 1 , 1 , 0 , 0 , 0 , 0 , 0 , 0 ), 
E7± = ±(-1, 0,-1, 0,0, 0,0,0), 
= ± 


1 1111111 
2 ’ ~ 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 ’ 2 


E± = ± 


( 1 

1 

1 1 

1111 

V 2 ’ 

~ Y ~ 

2’ 2’ 

2’ 2’ 2’ 2 


±( 0 ,- 1 ,- 1 , 0 , 0 , 0 , 0 , 0 ). 


(A. 8 ) 


In the main text we will need two particular linear combinations of these generators 


Qi — Hi ± 2 H 2 ± 2.^3 ± 2 Hi , Q 2 — H^ ± 2 Hi. (A.9) 


The other representations can also be identified. The roots in the representation (10, 5) 
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are the following ones 


/^5 — 
fJ-b ~ Oil = 




h5 

CXi 


a 2 — 0^3 — 


fJj^ — 0^1 — 0^2 — 0^3 — 0^4 — 


'11 1 

1 1 

1 

1 


2’2’ 2’ 

2’ 2’ 



2) ’ 

'1 11 

1 1 

1 

1 


2’ 2’2’ 

2’ 2’ 



2) ’ 

1,0, 0,1,1, 

0,0,0), 



'111 

1 1 

1 

1 


2’ 2’ 2’ 

2’ 2’ 



2) ’ 

' 1 1 

1 1 

1 

1 

1 1 

2’ 2’ 

2’ 2’ 



2’ 2 


(A.lOa) 

(A.lOb) 

(A.lOc) 

(A.lOd) 

(A.lOe) 


where we identified the hve 10 representations of SU{5)gut with their weight under 
SU{5)±. We called the highest weight of SU{5)± in the fundamental representation /is 
and the simple roots of SU (5)_l. We can apply the same procedure to the representation 
(5,10) and the result is 


/iio= (1,0, 0, -1,0,0, 0,0 ), 

'1 1 1 


hlO — Oi 2 — 

hio — Oil — a2 = 


11111 
2 ’ 2 ’ ~ 2 ’ ~ 2 ’ 2 ’ 2 ’ 2 ’ 2 
1 11 11111 
2 ’ ~ 2 ’ 2 ’ ~ 2 ’ 2 ’ 2 ’ 2 ’ 2 


hio — 0:2 — as — (0,1, 0, —1,0, 0, 0, 0 ), 
hio — ai — a2 — as = (0, 0,1, —1,0, 0, 0, 0), 


fiio — a2 — as — a4 — (0, 0, —1, —1, 0,0, 0, 0 ), 

fiiQ — Oil — 2a2 — as = 


111 11111 
'2’ 2’ 2’ ~2’ 2’ 2’ 2’ 2 

/iio - ai - a2 - as - a4 = (0, -1, 0, -1, 0,0, 0, 0), 

/iio — ai — 2 q; 2 — as — a4 = 


f 1 

1 

1 11111 

to 1 


2’ 2’2’2’2’2 


/iiQ — Oil — 202 — 2q;s — a4 — (—1, 0, 0, —1, 0,0, 0, 0) , 


(A.lla) 

(A.llb) 

(A.llc) 

(A.lld) 

(A.lle) 

(A.llf) 

(A.llg) 

(A.llh) 

(A.lli) 

(A.llj) 


where we called /iio the highest weight of the antisymmetric representation of SU{5)±. 
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B Details on Eg models 


In this appendix we gather some additional details regarding the models analysed in 
section We start by discussing how invariance under Es transformations constrains 
the possible couplings. We recall that the adjoint representation of Eg has the following 
branching rule when decomposing Eg under the maximal subgroup SU{5)gut x SU{5)± 

248 ^ (24,1) © (1, 24) © ((10, 5) © c.c.) © ((5,10) © c.c.). (B.l) 

It is convenient to introduce a basis of vectors ei,..., 65 for the fundamental representation 
of SU{5)± and this implies that a basis for the 10 representation of SU{5)± is given by 
Ci A Cj with i 7 ^ Since the 10 and the 5 representations of SU {5)gut sit in the 5 and 
10 of SU{5)± respectively we will label them using the basis vectors we just introduced as 
lOj and the Eij. Then invariance under Eg transformations boils down to invariance under 
SU{5)gut X SU{5)± transformations and therefore the following couplings are possible 

lOi ■ lOj ■ 5*^', (B.2) 

lOj ■ , ^ijklm 7 ^ 0 ) (®'2) 

where we have raised indices with the It is also important to analyse the possible 
couplings between the helds charged under SU{5)gut and the singlets that come from the 
adjoint representation of SU{5)±. Using again the basis of vectors e* of the fundamental 
representation of SU{5)± the elements in the adjoint representation will be e* © and 
labelling the helds in the (1,24) of SU{5) x SU{5)± as 1* the following renormalisable 
couplings with the helds charged under SU{5)gut are possible 


10 , 

■ 10^ ■ 1} , 

(BA) 

5p- 

' ■ li. 

(B.6) 


Now we will turn to a more detailed description of the Yukawa couplings present in the 
models analysed in section 

will denote the dual basis for the antifundamental representation of <S't7(5)_L as and similarly 
a basis for the 10 representation of SU{5)± is ef A et with i ^ j. 


50 



4+1 model 


In this model the Higgs held $ breaks SU{5)± to S{U{4:) x This induces the 

following breaking pattern for the fundamental and antisymmetric representations 


SU{5)^ - 

-A 5(f/(4) X H(l)) 

5 - 

—y 4 1 © 1 4 

5 5 

10 - 

—)■ 62 © 4 _3 


5 5 


(B. 6 ) 


where the subscript denotes the charge under the U{1) that comes from the trace of the 


f/(4) factor. Thus, given the unfolding (B.l), we have that the relevant sectors are those 
that appear in table A particular basis for these sectors is 


SU {5)gut 

S{U{4) X f/(l)) 

10 a 

4i 

5 

10 , 

1 4 

5 

5a 

62 

5 

5, 

4 3 


5 


Table 3: Different sectors for the 4+1 splitting. We do not show the conjugate represen¬ 
tations. 


10 sector 


5 sector 


10 „: 


5a ; < 


ei 

62 

63 

64 


61 A 62 
61 A 63 

61 A 64 

62 A 63 

62 A 64 

63 A 64 


, 


10 , 


> , 5, : < 


65 


61 A 65 

62 A 65 

63 A 65 

64 A 65 


(B.7) 


(B. 8 ) 
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It is immediate to see from the charge assignments that in order to have a SU{5)± 
invariant np-type Yukawa coupling we need to take the 10 m to be the 10a and the 5u as 
the 5a. Furthermore, the 5^ and 5 m need to be the 5a and 5;,, respectively. This comes 
from imposing a non-vanishing down-type Yukawa coupling as well as the possibility of 
having a correct chiral spectrum. Indeed, if we take the 5 m and 5u to be in the same 
curve, we cannot have three massless copies of 5 m and achieve doublet-triplet splitting. 
Note that using gauge invariance in the 5a sector we can set to zero all components of 
the wavefunction except for two which can be chosen to be the 62 A 63 and the 63 A 64 . 
In the following we shall the wavefunctions for these sectors and 02 respectively. The 
Yukawa matrices for the model are the following ones 


Yu = 


TT 


Yd/l — 


^PmP^ 


ttVi 

^PmP^ 


0 
0 

V 


0 


i Bx+bdy \ 

0 


/ 


0 

0 


0x-\-0y 


0 

‘^Pfi 

0 


0 01 + ^02^ y 

0x-\-b9y \ 

2 p^ 

0 


(B.9) 


3+2 models 


In this class of models the Higgs held $ breaks S'f/(5)j_ to S{U{3) x U{2)). In this case, 
the branching of the relevant representations reads 


^ S{U{3) X U{2)) 

5 ^ ( 3 , 1)2 © (3,1)_3 

10 ^ (3,1)4 ©(3,2)1 ©(1,1)_|, 


(B.IO) 


where the subscript denotes the charge under the trace of f/(3), which gives the helds in 
table m 

We pick the following basis 
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S'ly (5 )g[/t 

S(U(3) X U(2)) 

10 a 

(3,1)| 

10c 

(3,1)_| 

5 a 

(3.1). 

5c 

(3,2)1 

5 c 

(1,1)_| 


Table 4: Different sectors for the 3+2 splitting. We do not show the conjugate represen¬ 
tations. 


- 10 sector 


10 a: < 


ei 

62 

63 





(B.ll) 


- 5 sector 


5a : 


61 A 62 

62 A 63 > , 

6l A 63 

v / 


5f,: 


61 A 64 

62 A 64 

63 A 64 

< , 

61 A 65 

62 A 65 

63 A 65 


5 c : {64 A 65} . 


(B.12) 


In this class of models there are two possible choices for the 10m, namely either 10a 
or lOfc. Once the 10 m is chosen the 5f/ is uniquely hxed by demanding the possibility 
of having an up-type Yukawa. Then, for each of them, we have the freedom to choose 
between the 5 m and 5 d- The different possibilities are 


• 10m = 10a, 

• 10m = 10a, 

• 10 m = lOfc, 

• 10 m = lOfe, 


5c/ = 5 a, 
5c/ = 5 a, 
5c/ = 5 c, 
5c/ = 5 c, 


5m = 5f,, 
5m = 5c, 
5m = 5a, 
5 m = 5f,, 


5d = 5c. 


5d = 5c. 


5n = 5c. 


5 d = 5 a. 
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The Yukawa matrices for both cases in which we assign the IQm to the 10;, are the 
following ones 


Yu = 


TT 




0 

0 f 

Ox~\~^Qy 

/ 0 0 


0 

„ Qx~\~i>0y 
^pfl 

0 


. dx+b9y ^ 

0 


Yd/l — 


TT 


PraPri 


\ 


0 


0 2 e 


0 

{Gx )P/i 

3/2 

pm 




0 

- [(2 + 6)0xH“^^y]p^i 
■ 3/2 

Pm 

( Gb'^d^pl _ Wbdypl 


3p4 

i^n 




/ 


(B.13) 


For the cases in which we assign the 10 m to the 10^ the Yukawa matrix for the up quarks 


vanishes at tree level like it was already noticed in 23 , therefore we will not compute the 
Yukawa matrices for these two models. 


2+2+1 models 


In this class of models the Higgs held $ breaks SU{5)± to S{U{2) x f/(2) x Again, 


the representations of SU{5)± decompose as 


5H(5)x ^ S{U{2) X f/(2) X f/(l)) 


(B.14) 


5 ^ (2,1)>_. 

^ ^ 5 ’ 5 

10 -^ ( 1 , 1)6 _4 

^ ^ 5 ’ 5 


(1, 2)_2 3 © (1, 1)_2 _2 

^ ^ 5’5 ^ ^ 5’ 5 

(1,1)_4 6 © (2,l)l _4 © (1,2)_4 1 

^ ^ 5’5 ^ ^ 5 ’ 5 ^ ^ 5’5 


( 2 , 2 ) 11 
^ ^ 5’5 


Here the subscripts denote the charges under the traces of the two U{2) factors. The 
different sectors are displayed in table 

In terms of a particular basis we have the following 


- 10 sector 



10c : {^5} 


(B,16) 
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SU (5)g(7t 

S{U{2) X U{2) X U{1)) 

10a 


lOb 

(1,2)-|,| 

10c 

(1, !)-§,-§ 

5a 


5b 


5c 


5d 


5e 

(2.2)y 


Table 5: Different sectors for the 2+2+1 splitting. We do not show the conjugate repre¬ 
sentations. 


5 sector 


5 a : {ei A 62} , 5 b : {cs A 64} , 5 c : 


Cl A 65 
62 A 65 


5 . : 


63 A 65 

64 A 65 


, - < 


Cl A 63 
62 A 63 
Cl A 64 
62 A 64 


(B.16) 


In order to obtain the different models we proceed as before, namely we pick a par¬ 
ticular sector to be the 10 m and check whether one can choose a 5u such that a SU{5)± 
invariant coupling IOmIOmSc/ exists. Then, we check the different assignments of 5 m 
and that allow for a coupling 10 m5m5d- The different possibilities are 


10m = 10a, 5(7 = 5a, 5 m = Sfe, 5 d = 5c 


10m = 10a, 5(7 = 5a, 5 m = 5c, 5 d = 56 


10m = 10a, 5(7 = 5a, 5 m = 5^, 5d = 5e 


10m = 10a, 5(7 = 5a, 5 m = 5e, 5^ = 5^. 
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One can also consider the possibilities that arise from taking the ones above and permnting 
the two U{2) factors. However, these models are physically eqnivalent. Note that again 
in the 5e sector after using gauge invariance two components cannot be set to zero, and 
we will choose them to coincide with the roots ei A 64 and 62 A 64 calling them 0i and <p 2 
respectively. For the first model we find the following Yukawa matrices 


Yu = 


TT 


‘2pmP^ 


Yd/l — 


( 0 
0 

6x-\-h0y 


TT 


0 e 

, 0x~\~b0y 

0 


9x~\~b0y \ 

2 p/i 

0 


/ 


/ 


+ P^2) 


\ 


0 0 e 1/13 

0 e 1/22 e 1/23 
\ ^ 1/31 £ 2/32 1 + £ 2/33 


where 


2/13 = 

2/22 = 

2/23 = 

2/31 = 


9y {2bp^^ Y cp^^) Y (2pni Y ^^ 2 ) 

( 2 pni+Pi, 2 )^ 

9y (4 + 1) + 4(6c + l)py 2 Py^ + (P + 1) 

( 2 Pni+Pit 2 )^ 

3(6 — c)py_^Pfj_2 {9y {2bpyi^ + cpy,^) + (2piti + P/12) ^x) 

(2p/ii + P/ 12 ) 

6'y (26p^^ + cp^J (4 (6^ + 2) p^^ + 4(6c + 2)py^py,^ + (P + 2) pf^^) 


(2p/ii + P /12 


P32 — 


(B.17) 


2/33 — 


(26p^j^ + ^PfJ-P 
2p/ii Y P/12 

P^i^x (4 (6^ + 1) p^^ + 4(6c + 1)P/12P/11 + (c^ + 1) p^ 2 ) I 

( 2 p/ii + P/ 12 ) 

pIA ((3c - 26)p^ 2 + 26p/iJ (4 (6^ + 1) p^^ + 4(6c + l)p/i 2 P/ii + (c^ + 1) P^ 2 ) 

(2p/ii + P/ 12 ) 3pm^ 

6(fe - c)VmiPL (26p/ii + CP/. 2 ) + (2p/ii + P/ 12 ) Ox) 

(2p/ii + P/ 12 ) Vm 


(B.18) 


The result for the second model is presented in the main text in 


section For the third 
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model we find the following Yukawa matrices 


Yu = 


TT 


‘^PmP^ 


0 
0 

. 0x-\-i>0y 

\ 

/ 


r. dx+b8y \ 

2 p^ ' 

8x~\~b0y 


^ ^ 0 


0 


Yd/l = 


TT 


pm pa 


\ 


(B.19) 


0 0 0 
0 0 e 1/23 

\0 e 1/32 02 + e 1/33 J 


where 


1/23 — 


1/32 — 


(& - {(t)l0xPrh + 02Pp 2 (c^'y “ O^)) - &02PpP^^ (^^'y + “iOx 


3/2 3/2 
Pm Pm 


{h - c)p^J? (02Pm 2 (3c6'y + Qx) - c(l)i9yPfa) + c(j)2p^^p%^ {bOy - 9, 


■il2 


3/2 3/2 
Pm Pm 


(& - c)p^l^pfa (2hpl9xp%^ - p%pm^ {9x - 3c0y)) - cpypy,^p^jap^J^ (0^ - b9. 


Vss — — 0 r 


p^ p^- 
rmrm 


2cpy^9y (b{c - b)plp^J^ + {b + c)py^pyp)P + 3(6 - c)p2^p: 


2 , 3/2 

m 


Pm Pm 


2 pp 0 a; {^b‘^plp%^ + {b-c){b + c)py,^pyp^ja - cpiyja^ 


3 3/2 

PmPrh 


(B.20) 


Finally in the fourth model we face the problem that the matter curve where the sector 
5e is localised is actually a singular variety and therefore it is not clear how to correctly 
identify the various families for the 5 m sector. Since this constitutes a major issue in the 
analysis of the model we will not discuss it further. 


C Local chirality 

One important consequence of the addition of fluxes as already remarked in the previous 
section is the generation of chiral matter in 4d. In fact, while a non-zero vev for the Higgs 
background $ gives chiral matter localised in some matter curves S, this chiral matter 
lives in the six dimensional space x E and therefore the resulting four dimensional 
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spectrum will be non-chiral after dimensional reduction. Non zero fluxes threading the 
matter curves change the situation and if for a particular matter curve Sp 


Tr(F) ^ 0 


(C.l) 


one of the two chiralities for the matter living in Sp will remain in the massless spectrum 


in 4d. In our local analysis we cannot evaluate whether the condition (C.l) is met for a 
particular sector of the theory for this would require knowledge of Sp and the fluxes F 
globally in Sgut- Nevertheless we will be able to analyse the chiral spectrum of our model 


using the notion of local chirality. Local chirality was discussed in 25 -27 and it amounts 
to asking for localisation of the matter wavefunctions in a certain sector p around the 
Yukawa point. To understand concretely whether matter in particular sector is localised 
or not it is useful to look at the same problem in the T-dual setup of magnetised D9- 


branes. As explained in 24 the gauge field in the z direction is identified with and 
so Fxz = F>x^ and Fyz = Dy^. In the case of T-branes we additionally have the flux 
Fzz = <h^]. In this situation local chirality at the Yukawa point can be studied looking 

at the Dirac index around this particular point. For a representation TZ the Dirac index 

is 


indexT^T^ = 


48(27r) 


Tr 7 ^ F AF AF - -Tiv. F A TtR A R 


8 


7 ^ 


(C.2) 


and asking for local chirality at a particular point amounts to asking for a non vanishing 
integrand at this particular point. Since we take zero curvature in our model the quantity 
we need to evaluate is the following one 


I'll — {F A F A F)^-y-^- — i TTTi(^Fxx{Fyy, Fzz} + Fxz{Fyx-, Fzy} + 


(C.3) 


Fxy}Fyzi Fzx} }FxX-i Fyz} Fzy }Fxyy Fyx} Fzz {Fxz-, Fyy} Fzx} ■ 


As shown in [^, we hnd that the wavefunction in a representation TZ is localised when 
X-ji < 0. This can be evaluated for the the different sectors that appear in the Yukawa 
couplings which yields 


^IOm 

- 2m^c^qR{lQM) 

(C.4a) 

= 

- 2m^c^qR{%M) 

(C.4b) 

I^u = 


(C.4c) 

Ilo = 


(C.4d) 
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where we took 6=1 and /i m, as in section Using the charges in table [T] we hnd 
that we need to restrict the fluxes to 

Ml < 0, 4Ni < -Ny, 4N2 > -Ny (C.5) 

together with either 

3Mi < 2Ny < 0, 4M2 + Ny + 2 Mi > 0 (C.6 ) 

or 

0 < Ny ^ —-Ml, QM 2 — Ny 3Mi > 0. 


D Holomorphic Yukawa couplings computation 


In this appendix we present the explicit computation of the Yukawas via residues as 


explained in subsection 4^ Taking into account the non-perturbative corrections, the 
holomorphic Yukawa coupling is 


Y = ml'K'^fabc Resp [rj’^ifh^y] , 


(D.l) 


where 


r] = -*4) ^ [h^y + iedJody (<h ^h^y) - iedyOod^ ($ ^Ky)] • (D.2) 

Here is a different function for every sector. In particular, 4) is the action of (4)) on each 


sector, as in eq.(4.7), and hxy are given in (4.32). Finally, 9q is a holomorphic function on 


S'gut which we take linear as it appears in (3.3). Given this information we can proceed 


to compute the different r] functions which read. 
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10 m 


hio. 


= bx — y 


<o„/7 


Ml iWm 



1 m \ 

det$io^ 

V / 


+ 0(6^) 


10m 


w?x 


2fi^{9x + b9y)\ + 3 _^^ 3 _^ / 2mA 

(det$ioM)^ 

{9x + b9y) 3 _i, 3 _i ( 2mA(6-i) 


(D.3) 

(D.4) 




^5m = - ^o) -{y- yo) 

iri% HI = 

•OM' 'Om 


ml Hi * 

* 5m 


m 


(( 2 a + l)A + 2 /s:) - m^x \ {2a + l)X + 2K 


+ 6 


ml{b9y + 9x)}i 


2-i 

5m 


(2a + l)^^m (fcj + 2(3 - i)(bx + y)) 


(D.5) 

(D, 6 ) 


(detts.,)' \ (2K + ( 2 o+l)A)( 2 a+l)^=[fc 5 „ + (3-!)(6a + 9 )] 


+ 6- 


+ 6 - 


ml{b9y + 9x)hH / 4mm 

(det$ 5 ^)^ y (3 — i)m^x + [2 k, + (2a + 1)A](3 — i)K 
2ml *(2a + V)y,‘^{b9y + 9^) [(2a + 1)A + 2,k\ ( 2m [(2a + 1)A + 2k] 


(det4>5^ 


(det$ 


) 


2ml-H\2 / -2m((2a + l)A + 2K)m30j, 
(det<l>5^)3 1 -(det<l>5^)m^6'j; 


+ 0{^). 


(D.7) 


»[/ 


h^u/l^u — 1 


H^u/Hu = 


2X‘{9x + b9y) 


$ 


5[/ 


$3 

5(7 


+ 0{e‘) 


(D. 8 ) 

(D,9) 
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Hohln = 1 
= - 


1 2a^^{6x + b6y) 




— e 


$3 

5d 


+ 0{e^). 


(D.IO) 

(D.ll) 


Once we have these we can simply apply (D.l) to hnd the Yukawas. The computation 


of multivariate residues can be quite involved, however, in our case the matter curves 
intersect transversely, which means that these can be calculated in a straightforward way. 
Let us illustrate how this works. Consider the following residue at a point p G C^, 

f{x,y) 


R = ReSr: 


(D.12) 


where /, cxi and ct 2 are holomorphic functions. If ai and ct 2 meet transversely at p, 


meaning that the determinant of the Jacobian at p, 

d{ai,a2) 


Jio'l, 0'2)\p — 


d{x,y) 


(D.13) 


is non-zero, we can perform the following change of variables {u,v) = {ai{x,y),a 2 {x,y)). 
Then, the residue is 


R = Res,- 


</(o-i,o-2)r 


uv 




(D.14) 


where f{u,v) = f{x,y). Using this expression in our case leads to the result (4.33) in the 
main text. 


E Zero mode wavefunctions in real gauge 

In this appendix we discuss the computation of the real wavefunctions, both perturba¬ 
tive and non-perturbative. Although the structure of the couplings can be computed at 
the holomorphic level via residues, the normalisation factors and kinetic mixing neces¬ 
sarily depend on solving the D-term equation. Thus, as explained in the main text, we 
mainly use these solutions to compute the kinetic terms which is then combined with the 
holomorphic computation to give the physical Yukawa couplings. 

There are two different kinds of wavefunctions in our model, those that correspond to 
sectors which are charged under an Abelian subgroup of S{U{2) x U{2) x t/(l))_L and the 
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ones that transform as doublets of the first U{2) factor. These correspond to the Higgses 


and matter, respectively. The former have been computed in 24-26 and the latter in 26 


so we just sketch the computation and refer the reader to those for further details. 


Perturbative zero modes 

Higgs sectors 

These sectors are not charged under the T-brane background as can be seen from table 
Thus, we can compute their wavefunctions by applying the techniques in 24 26 . The 


zero mode equations (5.1) can be recast into a Dirac-type equation, namely 

\ / 


/ 


0 

Da, 

Dy 

D, 

-Da, 

0 

-D, 

Dy 

— Dy 

D, 

0 

-Dot 

-D, 

— Dy 

Ds 

0 


/ 


0 






u 


= 0 


(E.1) 


/ 


with 


= da, + -{qrX - qsy) Dy = dy- -{quy + qsx) = 2iy?{x - y) (E.2) 

and Dfn their conjugates. Here we used the following particular gauge 

A = -Qn^ydy — ydy — xdx + xdx) + -Qs{xdy — ydx + ydx — xdy) — -rri^c^Pi{xdx — xdx). 

(E.3) 

which reproduces the total flux, 

F = iQR{dy A dy — dx A dx) -|- iQs{dx A dy + dy A dx) + im^c^Pidx A dx. (E-4) 

The quantities q^i and qs are the constant flux densities that are shown in table and 
for concreteness we have taken the case of the up-type Higgs. The down-type Higgs can 
be obtained from this one by some simple replacements, see below. 


Following 24,25 one can solve this system of equations which yields 
/ a Al \ 

'^u = Iv Xu, 


\ 


_ Q^{\A'^-\y?)-is{^y+yx)+(x-y){Qux-Pu-Cu)y)) 




{P.h) 
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with Xu the lowest solution to 


~ + {QrY + {QsY)^u + — 0 


(E.6) 


and (u = shown in 26 we can multiply xu by an arbitrary holomorphic 


function of a particular linear combination of x and y and still satisfy the equations of 
motion. This function will only be determined once we impose boundary conditions that 
arise when we embed this model in a compact setup. However, the dominant contribution 
to the Yukawa couplings comes from the average value of such function around the Yukawa 
point, so we may approximate it by a constant. 


Now, looking at (4.7), we see that if we set k to zero, we can hnd the solution for the 


down-type Higgs by simply performing the following replacement on the solution above 


—y — djj,^. 


(E.7) 


Furthermore, the effect of k is simply a translation in the {x, y) plane which can be taken 
into account by shifting the solution 


x^ X-Xo, y ^y-yo. 


(E.8) 


This yields the wavefunction for Hd in a. gauge in which the vector potential is shifted 
with respect to the one used to solve the Hu sector. Thus, we need to perform a change 
of gauge for the background flux, namely. 


A{x - xo,y - yo) ^ A(x, y) + di/} 


(E.9) 


with 

V' = -QR{yoy-yoy-XQX+xox)-\--Qs{xQy-yox+yQX-xoy)--rn^c^Pi{xQX-XQx). (E.IO) 
Taking all this into account we hnd that. 


D — 


Id 


( j \ 

^ 2diJ? 

• Cd-^d 


V 


e - Xo, 2 /- 2 /o) 




with 


(E.ll) 


Xd ( x , v ) = e » 


— p^i\A‘^-\y?)-<is(^y+y^)+{^-y){C,DX-{\D-CD)y) 


(E.12) 
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Finally, is the lowest solution to 


- {M^ir + qR + ql)XD + 8d^iJ,^qs = 0 


(E.13) 


and (d = 


^pj^p-gR-qs) 

2{\p-qs) 


Matter sectors 

Unlike the previous sectors, these are charged under the T-brane since they transform 
as doublets of U{2). This makes the computation of the real wavefunction more compli¬ 


cated. In particular, as shown in 26 , it is not possible to find a simple solution in the 


case in which /x is non-zero since that leads to a pair of coupled partial differential equa¬ 
tions. However, in the limit /i ^ m these two equations decouple and can be effectively 
reduced to ordinary differential equations. In the following we present some details on 


the computation while a more detailed discussion can be found in appendix A of 26 


Similarly to what we did with the wavefunctions of the Higgses, we start by considering 
the 10m sector and we will obtain the 5 m by simple replacements. Furthermore, for 


simplicity we set 6 = 1 in the Higgs background (4.3). 


Since this sector is a doublet under the first U (2) factor, we need to have 

/ . \ 


tty 


(E.14) 


y ^xy J 


for which one can write an analogous equation to (E.l) where 


a = 



ip = 



(E.15) 


Then, following 25,26 , we first solve the F-terms equations (5.1a) and (5.1b) to write a 


in terms of </9, and then impose the D-term equation (5.1c) to hnd an equation for (p. 


Let us start by considering the case in which the primitive fluxes {Ep) are zero. Then 
solution to the F-terms reads 


= (h - 


ip = e 


{h - tn) 


(E.16a) 

(E.16b) 
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where ^ and h are doublets with components and and 


Pi = 


1 0 


^ = 


-IJ?{x — y) m 

-y^{x-y) 


ni?x 


From (E.16b) we obtain 


^ ^ (e — h) 


which can be used to write down the D-term equation for the fluctuations (5.1c) 


dxdsi + dydyi + dxfPdxi - {h - i^>C) = 0 


where we used that / does not depend on (?/, y) and we have dehned 


(E.17) 


(E.18) 


(E.19) 


m?xe 


A = = 

Finally, we can make the following change of variables 


IJ?{x — y) 

-y^{x-y) 


(E.20) 


U = 




e = {U - h) 


(E.21) 


and express (E.19) in terms of U 


dxdxU + dydyU - {dx^)'^~^dxU + dyU + dxf^P^~^dxU - = 0. (E.22) 

As advanced earlier, this is system of coupled partial differential equations for and U~ 


that have no simple solution. However, following 26 , these decouple in the limit m 3> /i. 


Furthermore, it is possible to show that there is no localised solution for U^, so we set it 
to zero. Then, near the Yukawa point pup = {x = y = 0} if we approximate / = logc + 
c^m?xx we hnd U~ = exp{Xxx)h where A the lowest solution to c^A^ + 4c^m^A^ —= 0. 
Taking this into account one hnds 


^ 10 + = xio 


/ ^ \ 
0 

V 0 / 




^ 10 - = 7^0 


/o\ 

0 

yi/ 


e 


(E.23) 


where = ^^m'^c^xxl 2 _ ^\xx with gj holomorphic functions of y. 
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Switching on the primitive gauge fluxes amounts to replacing dx,y —)■ in the D- 


term, where Dx,y is deflned in (E.2), and analogously for d in the F-terms. Then, taking 


/r m one finds a localised solution for U and the wavefunction reads 

/o^ 


/ ^Aio ^ 
2 


= 




\ 


_ »-^Cio 

m? 

0 




= 7lo 




0 

W 


e 


(E.24) 


where Aio is the lowest (negative) solution to 


m‘^(Aio - qn) + Ac^ {c^m^{qR - Aio) - A^g + + g|) =0 

and Cio = ~?s/(-^io ~ The scalar wavefunctions Xio are 


-\y\7-qs(=^y+ys)+Xiox(x-(^oy) g.^y ^ 


(E.25) 


(E.26) 


where gj holomorphic functions of y + Ciox, and j = 1,2,3 label the different zero mode 


families. Similarly to 13 we choose such family functions as 


9j = ml ^{y + Ciox) 


3-i 


(E.27) 


Now that we have the solution to the 10 m sector, we can obtain the corresponding 
solution for the 5 m by some replacements. Similarly to the case of the down-type Higgs, 


we arrive at the solution by performing the shift (|E.8|) where now 

Xq —)■ 0, 


V 

2/0 -, 
a 


(E.28) 


since we need to take the limit /i <C m as well as k m, keeping k/= v finite. Again, 
we have to perform a change of gauge with parameter 


~ % % 

"0 = ^Qnii^y/a - vy/a) + -Qs{,yxla - vxja) 


(E.29) 


which is the same as "0 but using the fact that, when taking fi,K m with = u 
finite, we have xq = 0 and yo = j a. Then, the wavefunction for the matter 5 m sector 
reads 

( As \ / 0 ^ 




75 


-iAAs 

m ^ 


V 


0 


/ 


0 

w 


e^^-f/Y,{x,y-v/a)E^ (E.30) 
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where A 5 is the lowest solution to 


"i^(A 5 - qn) + \^c^{c^m^{qR - Aio) - A 5 + + g|) = 0, (E.31) 

and C 5 = —qs /{^5 — qn)- Finally the wavefunctions xl are 

xi(x, y) = e^(l"l^-ls'l")-'?s(-W)+A5x(x-C55)^i ^ (E_ 32 ) 

where g\ are in ?/ + and i = 1,2,3 is a generation index. Analogously, the family 
functions are 

9l{y + C 5 X) = ml-\y + C,5xf~" ■ (E.33) 

Non-perturbative zero modes 

Here we discuss the correction to the real wavefunctions due to the non-perturbative 
effects. We start by considering the Higgs sectors and then move on to the matter wave- 
functions. 


Higgs sectors 

The analysis of the correction to these two sectors are identical to those that appeared in 


section 5.1 of p5) and in appendix A of 26 . The zero mode equations read 


d{A)a 

= 0 

(E.34) 

a(A}(^ + i[($),a] +e^eo^^l^A)a 

= 0 

(E.35) 

uAd(^A)a - 

= 0 

(E.36) 


that reduce to the following in holomorphic gauge for the up-type Higgs 


dstty - dytts = 0 (E.37) 

^fh^xy y'jOjfh, If^y^x^rh ^x^y^fh\ T ) (E.38) 

regarding the F-terms while the D-term reads 

{dx + xqR - yqs] + {dy - yqR - xqs] ay 2i^^{x - y)ip^y (E.39) 

= ie9^ {yqR xqs} (fxy - ieOy {-qRX + yqs} (fxy 
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The first order correction to the wavefunction is then 


= rn,^^e^^-y)iQus-{Xu-i:u)y))-[^ 


(E.40) 


with A[/, Cf/ dehned as in (E.5) and 

= ■^{Cux-{Xu-Cu)yf{Ox + Oy) + ^-^{x-yf + ^{x-y){Cuy + {^u-Cu)x) (E.41) 

where the constants and 62 are 

2/^2 


= 

^2 = 


X2 


{(^x{(1r{Cu — Xu) + qsCu) + (^yi^nCu — QsiCu — Xu))} (E.42) 


2/i^Ct/ 

X2 


{dxiqn + qs) + Oy{qR - qs)}- 


(E.43) 


Notice that the holomorphic terms in Tu, those that depend on and 62 , are there to 
satisfy the corrected D-term equation. 

Similarly to the tree-level case, the corrected wavefunction for the down-type sector 


can be obtained from this one by performing the replacements (E.7 E.8) and taking into 


account the change of gauge (E.9) so we do not write the explicit result. 


Matter sectors 

Recall that the wavefunctions in the real gauge are used essentially to compute the nor¬ 
malisation factors and kinetic mixing since the structure of the Yukawa couplings can 
be computed via residues. For the matter sectors in our model, one can prove that the 
mixing and normalisation factors remain unchanged when including the non-perturbative 
effects just by analysing the structure of the equations. In the following we discuss this 
point and obtain the structure of such correction for the matter sectors. 

Just as in the perturbative case, let us start by turning off the primitive fluxes. Also, 
we consider the 10 m sector and obtain the 5 m by performing replacements. First, we see 
that the F-terms in the real gauge are solved by 

O' = ydi (E.44a) 

if = g{h — — edOo Xd^) = gU dx A dy (E.44b) 

with 

e //2 0 \ 

0 e -//2 ) 
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and <l> is given by $|iom (4-7), where we dropped the subscript 10 m for notational 


convenience. The doublet U can be expanded in e 


U = + ef/W + 0{e^) 

where in the solution for e = 0 that appeared in the last section, namely 


(E.46) 


jjio) ^ 


f/i°^ = 0. 


(E.47) 


Then, one may solve for ^ from (E.44b) as 


^ = ^( 0 ) + [^( 1 ) + dJody^^^ - dyOod^^^^^ + 0{e^) 

^(0) _ i<|.-i(f/(o) _/i) 


(E.48) 


and then solve for by plugging in this expression into the D-term for the fluctuations 


(|5.1c|). This yields = 0, in the limit y <^m. Thus, we find the following structure 


e+ = + 0 + O(e^) 


e- = 0 + ed'^ + O(e2). 


(E.49) 


Since determines both a and (p through (E.44), we find that this particular structure 


implies that the solution (E.24) looks like 

/o\ 

+ e 


"^10+ ~ 




0 

V/ 


+ 0{e^) ^ 10 ^ = 


/o\ 

0 

v*y 




w 


+ 0{e^) (E.50) 


Finally, it is possible to show that this structure still holds after including the non¬ 
primitive fluxes. 
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